w 


1 


& 


Atmospheric  Transmittance  Study  with 

“-a:  fiaatsr «- 


P art 

jnWWN  TO  Ticaj  niGIT/lL 

HAT  A niTTPThATTnil  rww  rmm 


SOaTinN  70  THE  jgg,  JTIGITAL  , 
JflTA  REGISTRATION  PROBLEM  ♦ 

'-—rrz. ^ / 


D D Ov 

inMSEIMI2|n 

MAY  2919^  U 

^ikgsEQirEUy 


% 


Contract 


— if 

-7  FiNAL/ba^  \ TJ/1/7F  (=tfm>,-7f-rr-y-.  PT-3  / 

i <$J  1 s-  3$ 


by 


Sandra  K./ Weaver 
Rufus  E.y Bruce 

Joseph  H. jpierluissi 


Prepared  for: 


United  States  Army  Electronics  Command 
Atmospheric  Sciences  Laboratory 
White  Sands  Missile  Range 
New  Mexico 


Submitted  by 


Physics  Department 
and 

Electrical  Engineering  Department 
The  University  of  Texas  at  El  Paso 
El  Paso,  Texas 

yot  r?? 


Rufus  E.  Bruce 
Joseph  H.  Pierluissi 
Project  Directors 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DDC  CONTAINED  A SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


REPORT  DOCUMENTATION  PAGE 


SECURITY  CLASSIFICATION  OF  THIS  FAOE  (When  Dote  entered) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  GOVT  ACCESSION  NO.  S.  RECIPIENT'S  CATALOG  NUMBER 

PR4-76-DC-30  * 


REPORT  NUMB 


4.  TITLE  (end  Subtitle) 


Solution  to  the  SMS  Digital  Data 
Registration  Problem 


7.  AUTHORfaJ 

Sandra  K.  Weaver 
Rufus  E.  Bruce 
Joseph  H.  Pierluissi 


#.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Physics  and  Electrical  Engineering  Depts. 
The  University  of  Texas  at  El  Paso 


S.  PERFORMING  ORG.  REPORT  NUMBER 


S.  CONTRACT  OR  GRANT  NUMBERS 


DAEA18-76-C-0019 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  • WORK  UNIT  NUMBERS 


MHi/rrrn 


11.  CONTROLLING  OFFICE  NAME  AND  ADORESS  12.  REPORT  DATE 

United  States  Army  Electronics  Command  October  31,  1976 

Atmospheric  Sciences  Laboratory  <*•  number  of  pages 

White  Sands  Missile  Ranee,  New  Mexico  88002  82 


MONITORING  AGENCY  NAME  • ADDRESS (II  dl  Heron!  from  Controlling  Ollteo)  IS.  SECURITY  CLASS,  (ol  tbit  report) 

Unclassified 


ASSI  FI  C ATI  ON  / DOWN  OR  ADIN  G 


IS.  DISTRIBUTION  STATEMENT  (ol  thlt  Report) 


Approved  for  Public  Release:  Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (el  the  ebttrect  entered  In  Block  70,  II  dllte rent  treat  Report) 


IB.  SUPPLEMENTARY  NOTES 


Contracting  Officer  Technical 
Representative 


19.  KEY  WQftOf  (Contlnum  on  mwtm  ff*  II  and  identity  by  bit 


Richard  B.  Gomez 

Atmospheric  Sciences  Laboratory 

White  Sands  Missile  Range 


T'iTV  Tr.’.Lli 


Satellite  Data,  SMS  Data,  Automatic  Landmark  Scanning 


If  (Wtmwy  MB 4 Idamttfy  by  blaab 

The  method  of  treatment  as  well  as  the  method  of  solution  of  the 
registration  of  SMS  digital  data  are  delineated  herein.  This  portion 
of  the  final  report  concerns  task  3.2.7.  V 


Foreward 


This  is  Part  III  of  the  final  report  under  Contract  DAEA18-76-C-0019 
entitled  Atmospheric  Transmittance  Study  with  the  Meteorological  Satellite 
Technical  Area  of  the  Atmospheric  Sciences  laboratory  at  White  Sands  4^0'  ^ ^ 
Missile  Range . Part  I contains  a study  and  development  of  band  models 
for  use  in  connection  with  techniques  for  the  calculation  of  atmospheric 
transmittance  along  slant-paths.  Part  II  contains  the  report  on  the 
studies  related  to  the  inversion  of  the  radiative  transfer  equation  for 
temperature,  composition  and  possible  cloud  correction  techniques  in  the 
15y  COj  band  region.  Also  included  there  is  a discussion  of  the  method 
used  in  this  study  for  the  calculation  of  atmospheric  transmittances 
using  line  spectral  parameters.  Part  III  deals  with  the  study  of  SMS 
digital  data  and  their  use  in  severe  storm  and  cloud  studies. 
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1.  Discussion  of  Effort 

At  the  onset  of  the  contract,  discussion  of  the  objectives  brought  to  light 
the  necessity  of  an  accurate  SMS  digital  data  registration  scheme  to  complete 
task  3.3.7,  an  examination  of  the  relation  of  severity  of  storm  to  cloud 
development.  Since  an  interest  in  that  area  was  expressed  by  the  contractor 
and  it  had  been  stated  in  the  contract  that  the  problem  of  registration 
would  be  looked  into,  it  was  decided  to  tackle  the  problem.  Many  other 
organizations  have  already  and  are  dealing  with  the  problem.  It  would 
ultimately  be  our  task  to  assimilate  the  knowledge  and  expertise  of  those 
individuals  in  the  field  of  registration  to  arrive  at  an  accurate  scheme 
adapted  to  the  needs  of  the  contractor. 

A straightforward  geometrical  approach  to  the  problem  was  proposed 

by  the  contractor.  While  work  began  on  the  dynamical  aspects  of  geometrical 

transforms  and  Euler  angles,  a survey  of  the  field  was  conducted.  It  was 

discovered  that,  among  others,  E.  Smith,  now  of  Colorado  State  University 

(C.S.U.)  and  formerly  of  Wisconsin,  and  C.E.  Velez,  of  NASA  Goddard  Space 
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Flight  Center  (NASA  GSFC)  are  experts  in  the  field  of  registration  and 
that  their  schemes,  the  former  of  a short-term  attitude  predictive  nature 
and  the  latter  of  a long-term  orbit  and  attitude  predictive  nature,  arc 
representative  of  the  current  two  approaches  to  the  problem.  Contact 
with  the  two  was  initiated.  A trip  to  CSU  to  talk  to  E.  Smith  was  arranged. 

A demonstration  of  CSU's  registration  capabilities  and  a copy  of  a write-up 
on  the  system  E.  Smith  and  D.  Phillips  had  developed  at  Wisconsin^was 
provided  along  with  an  offer  of  a copy  of  the  software.  Further  study  of 
the  write-up  indicated  that  the  software  would  indeed  be  helpful  so  the 
request  was  made.  A copy  of  the  software  was  received  in  early  December. 
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Dr.  Velez  was  also  contacted  and  offered  to  send  us  a copy  of  a report 

on  the  system  his  group  at  NASA  along  with  a group  from  NOAA  National 

3 

Environmental  Satellite  Service  (NOAA  NESS)  had  developed.  We  requested 
a copy  of  the  software  but  this  was  not  seen  to  be  feasible  because  their 
package  was  actually  part  of  two  huge  multipurpose  software  packages 
that  require  a lot  of  storage  and  that  are  virtually  undocumented.  Dr.  Velez 
suggested  a trip  to  NASA  GSFC  sometime  after  the  first  of  the  year  to  determine 
if  their  system  would  be  suited  to  the  needs  of  the  contractor.  Also,  his  group 
was  in  the  process  of  rewriting  and  consolidating  their  software  into  one 
package  for  their  new  PDP  1170,  which  is  compatible  with  Met.  Sat's  PDP  1145, 
and  he  anticipated  the  completion  of  the  program  by  January.  A trip  to 
Washington,  D.  C.  to  include  visits  to  both  NASA  GSFC  and  NOAA  NESS  was 
arranged . 

The  trip  to  Washington,  D.  C.  on  January  5-6,  1976  provided  invaluable 
insights  into  the  registration  problem  (see  trip  report  in  AppendixA]. 

Based  on  the  facts  that  1)  the  adaptation  of  Velez's  scheme  to  the  PDP  1170 
was  not  completed,  nor  2)  was  it  accurate  enough  for  our  purposes  (30-40  pictel 
elements  vs.  the  needed  1-2  pictel  elements),  plus  the  fact  that,  3)  our  attempt 
at  a straightforward  geometrical  scheme  proved  to  be  futile  [See  Summary  of 
SMS  Data  Processing  System  in  Appendix®],  it  was  decided  all  around  that 
the  best  route  would  be  to  adapt  a version  of  Wisconsin's,  of  the  D.  Phillips- 
E.  Smith  routine,  a copy  of  which  E.  Smith  had  already  sent  us  and  the 
accuracy  of  which  was  professed  to  be  1-2  pictel  elements.  It  was  also 
decided  that  it  would  help  if  a copy  of  the  mini -computer  adapted  Wisconsin 
version  in  use  and  offered  by  J.  Billingsley's  group  at  NASA  GSFC  was 
obtained.  A letter  was  sent  requesting  the  software  and  documentation  of  the 
navigation  segment  of  their  METPAK.  Meanwhile,  an  attempt  was  made  to 
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decipher  NOAA  NESS's  versiojj,of  their  registration  program,  a copy  of 
which  was  received  during  the  trip  to  Washington.  Because  the  large 
program  is  almost  totally  undocumented,  it  was  decided  that  the  great 
amount  of  effort  that  would  be  required  to  comprehend  the  program  would 
not  be  justified  by  the  result.  Therefore,  we  proceeded  with  an  in-depth 
analysis  of  the  theory  behind  the  D.  Phillips-E.  Smith  routine  which 
necessitated  periodical  communications  with  E.  Smith.  A mathematical 
overview  of  that  analysis,  based  on  the  wtite-up  of  the  scheme  by  D.  Phillips 
and  E.  Smith,  follows  in  Appendix C The  requested  NASA  GSFC  software, 
which  is  a version  of  Wisconsin's  Me  IDAS  and  which  had  already  been  adapted 
to  a PDP  1145  system,  was  received  mid-April.  Based  on  the  fact  that  it 
had  already  been  adapted  to  a PDP  1145  and  that  its  documentation  was 
a great  improvement  over  that  in  E.  Smith's  program,  we  decided  to  use 
the  NASA  GSFC  version,  developed  by  C.  T.  Mottershead  of  the  Computer 
Sciences  Corporation  (CSC),  as  a basis  for  the  prototype  to  be  developed 
for  the  needs  of  the  contractor. 

Because  the  theory  behind  the  Mottershead  version  is  basically  that 
of  the  D.  Phillips-E.  Smith  scheme,  the  analysis  in  Appendix C also 
applies  to  Mottershead 's . A thorough  flowchart  and  investigation  of 
parameters  used  was  undertaken  in  the  process  of  which  correspondence  with 
Mottershead  was  and  is  still  being  maintained.  Basic  10  handling  routines, 
which  are  machine-dependent,  are  the  major  obstacle  to  the  utilization  of 
Mottershead's  scheme  as  is.  Approval  for  the  use  of  Met  Sat's  PDP  114!., 
in  conjunction  with  one  of  their  personnel,  on  which  to  adapt  the  navigation 
program  was  given  in  June  by  the  contractor.  However,  the  PDP  was  not 
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ready  for  routine  use  (thorough  hardware  and  software  system  check-out 
was  necessary)  until  recently.  In  the  interim,  sample  calculations  of 
various  segments  of  Mottershead's  navigation  program  (based  on  the  test 
landmark  data  furnished  within  the  program)  were  performed  by  hand  to 
determine  expected  values  of  parameters.  Time-consuming  iteration  calculations 
necessitated  computer  assistance.  Therefore,  the  process  of  revising 
Mottershead's  program  t i fit  UTEP's  IBM  capabilities  was  undertaken, 
i.e.  editing  capabilities  (mid-stream,  interactive)  were  removed  and  PDP 
system-based  software  (data-manipulative)  was  either  removed  or  revised. 

The  error-checking  process  is  still  on-going;  however,  an  undocumented 
listing  of  the  program  in  its  present,  unusable  state  is  included  for 
reference  [see  Appendix D ] . Mottershead  is  in  the  process  of  revising  and 
consolidating  NASA  GSFC's  navigation  program  for  use  on  their  new  PDP  1170. 

When  it  is  operational,  he  will  send  us  a card  deck  of  the  final  program 
to  be  implemented  on  Met  Sat's  PDP  1145. 

Further  consultation  with  Mottershead  will  be  necessary  in  order  for 
us  to  adapt  the  editing  capabilities  of  his  program  to  the  peripherals  of 
Met  Sat's  PDP  1145.  Also,  it  will  be  necessary  for  us  to  work  with  the 
engineers  of  the  Met  Satellite  Technical  Area  to  develop  the  needed  data 
10  package. 

Hie  organizations  with  which  we  had  been  in  contact  all  routinely 
use  video  refresh  capability  (CRT  screen  with  cursor)  to  locate  the 
landmarks  needed  to  determine  satellite  attitude  and  reference  satellite 
position  in  each  photograph  [see  Appendix  A].  Since  the  contractor  does 
not  have  the  necessary  video  refresh  hardware,  some  other  method  of  locating 
landmarks  needed  to  be  developed.  A group  at  NOAA  NESS  is  experimenting 
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with  a cross-correlation  scheme  [see  Appendix  A]  and  provided  us  with  a 
partial  copy  of  their  unfinished  software,  totally  undocumented.  Their 
scheme  works  best  with  visible  data  and  requires  at  least  a year  of  testing 
with  sample  data  because  it  is  based  on  having  predetermined  blocks  of 
expected  data  in  storage.  It  is  our  impression  that  a more  generalized 
scanning  method  applicable  to  both  IR  and  visible  data  would  best  suit 
the  needs  of  the  contractor. 

In  order  to  get  a better  idea  of  the  scanning  method  needed,  we  decided 
to  take  a look  at  a few  examples  of  dumps  of  digital  SMS  data.  We  looked 
first  at  IR  data  because  of  the  lesser  quantity  of  data  necessary  to  view 
a relatively  large  area  and  because  of  our  ultimate  aim  of  relating  IR- 
derived  temperatures  to  cloud  height  for  both  the  severe  storm  case  study 
and  incorporation  into  VTPR  retreivals.  By  looking  at  the  corresponding 
IR  laser  image,  we  were  able  to  define  approximate  count  boundaries 
between  land  and  ocean  or  cloud.  We  then  correlated  these  approximate 
count  ranges  for  land,  oceans  and  clouds  with  the  temperatures  represented 
using  the  chart  shown  in  Figure  1 and  making  sure  that  the  temperatures 
"made  sense"  with  what  we  expected  climatologically. 

Only  64  gray  shades  are  available  for  our  use  in  the  IR  because 
the  least  significant  bit  is  dropped  from  the  data  so  that  it  will  fit  ^ 

on  a 7 track  tape  to  make  it  compatible  with  the  UNIVAC  computer.  (This 
makes  it  necessary  to  multiply  all  count  values  by  4 to  compare  with  the 
chart  value  in  Figure  1.).  When  the  PDP  system  is  complete,  it  will  be 
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equipped  with  a 9-track  tape  drive  along  with  the  capability  of  providing 
all  256  gray  shades  in  the  IR. 


We  are  very  much  aware  of  the  fact  that  the  amount  of  land-ocean  contrast 


and  the  actual  boundary  counts  found  are  functions  of  the  landmark  location, 

the  time  of  day  and  the  time  of  year.  For  this  reason  we  proposed  a yearly 

study  of  land-ocean  contrasts  to  determine  which  landmarks  are  best 

for  different  times  of  the  day  and  year.  The  maximum  contrast  would  seem 

to  occur  during  late  summers  and  late  winter  and  during  mid-afternoon  and 

early  morning  since  at  these  times  are  maximum  land  temperature  extremes 

and  also  c^ean  surface  temperatures  change  relatively  little  diumally  or 
S' 

seasonally.  Because  at  certain  times  of  the  day  and  of  the  year,  there 
is  only  one  count  difference  between  land  and  ocean  it  is  sometimes  very 
difficult  to  determine  where  the  water  actually  stops  and  the  land  starts. 

(Itis  possible  that  when  all  256  gray  shades  are  made  available  to  us  a 
larger  relative  temperature  difference  between  land  and  water  might 
become  evident.).  Also,  because  the  resolution  of  the  1R  is  4 km  X 8 km,  each 
data  point  is  an  average  of  the  temperatures  in  that  block.  The  coastline 
could  actually  be  anywhere  within  that  block  of  data  so  the  uncertainty 
in  determining  a coastal  landmark  point  is  at  least  one  data  block, 
not  withstanding  any  calculatory  manipulations  of  data.  For  these  reasons, 
it  has  become  our  idea  that  the  best  way  to  approach  the  automatic  landmark 
retrieval  problem  is  with  a scheme  that  first  scans  in  the  IR  to  locate 
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a landmark  point  ♦_  one  pictel  element  and  then  to  call  up  the  corresponding 
visible  lines  and  elements  that  correspond  to  the  IR  line  and  element 
plus  uncertainty  (visible  data  resolution  is  1 km  X 1 km  so  each  IR  pictel 
element  is  equivalent  to  32  visible  pictel  elements)  and  scan  for  the 
landmark  in  the  corresponding  visible  block  of  data.  It  is  sometimes 
difficult  in  visible  digital  data  dumps  to  distinguish  between  the 
brightness  values  associated  with  low  clouds  and  those  associated  with 
highly  reflective  land  or  ocean  surfaces.  By  scanning  in  the  IR  first  where, 


except  in  high  latitudes  in  the  winter,  the  cloud  counts  will  most  certainly 
be  higher  (representing  colder  temperatures)  than  land  or  water  counts, 
we  have  the  capability  of  automatically  ruling  out  cloud-covered  landmarks. 

During  the  year  a "Landmark  Scanning  Program"  was  written  and  is  in 
the  course  of  testing.  The  program  is  designed  to  demonstrate  the  feasibility 
of  identifying  specific  coastal  landmarks  using  computer  techniques.  It 
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determines  the  coastal  outline  by  identifying  high  contrast  regions  with 
a first  order  difference  technique.  A cloud  identification  is  incorporated 
into  the  program. 

The  program,  included  in  Appendix  E,  has  not  been  optimized  and  should 
not  be  considered  as  a finished  product.  Several  cleanup  problems,  all  of 
a relatively  minor  nature,  must  be  completed.  An  example  of  these  problems 
is  the  fact  that  the  coastline  is  generally  displaced  eastward  and  southward; 
this  is  purely  the  result  of  the  differencing  technique  being  used. 

In  its  present  form,  the  program  demonstrates  that  coastal  features, 
including  islands,  can  be  efficiently  obtained  from  the  IR  data.  Although 
we  have  not,  as  yet,  obtained  the  desired  accuracy  of  +_  1 pictel  element, 
we  can  see  no  reason  why  this  will  not  be  accomplished  when  the  program 
cleanup  is  completed. 

The  present  technique  requires  the  computer  to  identify  the  coastline 
and  then  to  obtain  its  most  westward  point.  It  is  clear  that  another 
coastal  landmark  identification  will  require  a different  search  criteria. 

Landmark  identification  schemes  of  this  type  will  require  that  particular 
search  algorithms  must  be  associated  with  each  different  landmark.  This 
problem  may  be  overcome  by  obtaining  a more  generalized  pattern  recognition 
scheme  or  by  developing  a more  general  search  routine.  With  respect  to 


the  former  we  have  initiated  efforts  to  use  both  fourier  and  mellin 
transforms  in  the  recognition  scheme.  These  efforts  are  in  the  first 
stages  of  programming.  The  later  approach  will  not  be  pursued  until 
we  have  developed  successful  search  routines  for  several  different  types 
of  coastal  landmarks. 
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2.  Task  3.2,7 

The  contract  stated  in  task  3.2.7  that  a study  would  be  performed 
in  which  we  would  "examine  conventional  synoptic  data  and  satellite 
images  of  severe  storm  systems  and  determine  the  correlation  that  exists 
between  cloud  type  and  development  and  time  of  greatest  severity."  We 
were  unable  to  perform  this  study  because  of  unforeseen  problems. 

In  the  first  place,  in  a severe  storm  study,  the  nature  of  the 
image-making  process  necessitates  using  the  digital  data  from  which  the 
images  are  made  to  get  necessary  detail.  Conclusions  about  the  relative 
brightness  levels  of  clouds  in  the  visible  data  or  about  the  cloud  top 
temperatures  in  the  IR  cannot  be  made  from  satellite  images.  Two  factors  are 
the  basic  reason  for  this:  1)  the  film  density  is  not  constant  within 
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one  image,  let  alone  from  one  image  to  the  next,  and  2)  often  the  same 
exposure  setting  from  one  image  to  the  next  cannot  be  utilized  because 
of  such  factors  as  the  sun  angle  over  the  area  of  interest. 

Secondly,  in  order  to  use  the  digital  SMS  data,  it  must  be  registered 
accurately;  and,  as  has  been  shown  above,  this  is  no  easy  task.  Often 
landmarks  from  which  to  register  the  data  cannot  be  found  near  the  severe 
storm  cell  in  question,  especially  if  the  cell  is  enmeshed  in  the  cloud 
mass  of  a much  larger  system.  Because  1)  the  SMS  satellites  spin  as  they 
record  the  data,  each  pictcl  element  having  been  taken  at  a different  time, 
and  2)  the  satellites  are  not  perfectly  geostationary  and  their  positions 
must  be  mathematically  derived,  a simple  interpolation  between  the  pictel 
elements  of  known  landmarks  to  locate  severe  storm  cells  becomes  impossible. 

The  registration  must  be  accurate.  Overshooting  tops,  a sign  of  particularly 
active  and  well -developed  cumulonimbus  or  severe  storm  cells,  have  been 
seen  on  SMS  images  to  cover  only  a fraction  of  the  area  of  the  underlying 
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TRIP  REPORT 


I 


Kanal  of  the  University  of  Maryland,  who  has  worked  with  SMS  data  in  pattern 
recognition.  We  are  in  the  process  of  searching  the  literature  for  possible 
articles  by  Dr.  Kanal  on  that  subject. 

On  the  severe  storms  aspect,  NASA  GSfC  has  already  implemented  a 
scheme  in  their  METPAK  for  the  IDAMS  which  computes  divergence  and  vorticity 
using  interpolated  wind  fields  'from  SMS  digital  dataand  which  R.  Adler 
and  J.  Billingsley  also  demonstrated.  A vertical  velocity  scheme  is  planned. 

NASA  is  aiming  toward  a real-time  data  system  for  severe  storm  research.  J 

R.  Adler  and  C.  Peslen,  also  STORMSAT  researchers,  are  interested  in  , 

maintainig  close  ties  and  offered  their  assistance  when  needed.  R.  Adler 
is  also  interested  in  January  10,  1975  data  and  offered  to  track  down  B.T. 

Miers  request  if  the  data  had  not  yet  arrived.  It  had  not  arrived,  so 
R.  Adler  was  contacted  via  phone  January  13,  1976.  He  promised  to  look 
into  the  data  request  and  who  to  contact  at  NASA  about  being  placed  on 
the  orbit  parameter  mailing  list.  He  was  contacted  again  January  29,  1976 
and  said  that  about  half  of  the  January  10,  1975  tapes  were  on  their  way 
along  with  some  hard  copy  of  the  data  he  had  requested  on  his  own  be  sent 
to  Met  Sat.  Also,  he  suggested  we  ask  Dr.  Velez  about  the  orbit  parameter 
mailing  list.  At  NOAA,  R.  Gurd  demonstrated  the  capabilities  of  their  Man 
Machine  Interactive  Processing  System,  the  most  impressive  features  being 
the  VTPR  and  NMC-based,  large  scale  cloud  height  determination  schemes. 

M.  Young,  head  of  the  winds  section,  offered  some  advice  on  cloud  height 
determination  for  severe  storms  applications. 

Lastly,  the  quality  of  SMS  data  was  investigated.  J.  Lienesch  of 
NOAA  discussed  the  quality  of  IR  data  in  particular,  possible  limb  effects 

in  the  IR,  and  data  quality  control  done  by  NOAA.  He  also  discussed  his  % 

experience  in  working  with  IR  data  as  to  land~ocean  and  cloud  top  contrasts. 

He  suggested  H.  Jacobowitz  , also  of  NOAA,  be  contacted  about  limb  effects 
in  the  visible  data.  Dr.  Jacobowitz  said  that  such  a study  was  being  planned. 

As  to  the  orbit  parameters  sent  with  the  housekeeping  data.  Dr.  Velez  of 
NASA  said,  the  beta  values  are  relatively  accurate  and  require  only  minor 

! 

corrections.  J.  Ellickson  of  NOAA  said  they  will  be  working  with  NASA  on 
much  improved  and  many  more  orbital  parameters  being  sent  as  housekeeping 
data  by  the  end  of  '76.  He  provided  copies  of  two  SMSA  and  B data  manipulation 
reports . 
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Department  of  rhyrtet 
lit.  1'aso,  Tkxak  7!>!M1H 


Um'wrii^  of  fC*?**  **  £l  (P*6o 


March  2,  1976 


Commanding  Officer 
Atmospheric  Sciences  Laboratory  f 
White  Sands  Missile  Range 
White  Sands,  New  Mexico  88002 

ATTN:  Dr.  Richard  Gomez 

RE1:  Contract  DAEA  18-76-C-0019 

Dear  Dr.  Gomez: 

Attached  is  a short  summary  of  the  SMS  Data  processing  system  on  which 
Mrs.  Weaver,  Mr.  Guard  and  I are  working. 

This  report  outlines  our  objectives,  approach  and  to  some  degree  the 
flexibility  that  we  are  intending  to  put  into  the  system.  Before  a 
get  too  far  into  developing  this  concept,  I believe  that  it  should  be 
reviewed  by  you  and  those  personnel  at  the  Atmospheric  Sciences  Laboratory 
who  are  most  concerned  with  the  work. 

Would  you  please  advise  me  if  our  approach  is  satisfactory  for  your  needs. 
Very  truly  yours. 


Rufus  E.  Bruce 
REB/gla 


SMS  DATA  PROCESSING  SYSTEM 


Overall  Program  Concept 

A.  Input  Consideration 

B.  Projected  Output  Capabilities 

C.  Computational  Technique 

Current  Areas  of  Activity 

A.  Pattern  Recognition 

1.  Initial  programing,  goals 

2.  Results  of  current  testing 

3.  Further  desirable  capabilities 

B.  The  Registration  Transformation 
1»  Direct  geometric  approach 

a.  initial' concept 

b.  diagram  and  equations 

c.  result  of  testing 

2.  Eric  Smith's  registration  program 
Research  In  the  Near  Future 

A*  Modification  & Incorporation  of  Eric  Smith's  Work 

B.  Final  Landmark  Acceptance  Criteria 

C.  Severe  Storms  Applications 
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PROGRAM  CONCEPT 

A system  is  proposed  which  processes  SMS-GOES  visible  and  Infra-Red  data 
to  be  used  in  severe  storm  case  studies  and  atmospheric  radiation  transmission 
research.  Our  support  toward  the  overall  system  includes  development  of  soft- 
ware for  landmark  registration,  and  subsequent  registration  and  transformation 
of  desired  data  blocks  based  on  known  and  calculated  parameters.  Presently 

I * 

efforts  have  emphasized  techniques  for  handling  I.R.  data. 

This  software  is  currently  being  written  for  and  tested  on  the  UNIVAC  1108 
computer  at  WSMR.  Consideration  is  being  made  during  coding  and  documentation 
to  allow  implementation  of  the  final  system  on  a PDP-11  with  minimal  difficulty. 

Landmark  registration  will  be  automatic,  with  manual  decision  override 

I 

capabilities  during  critical  stages  of  processing.  Once  a suitable  number  of 
landmarks  has  been  identified,  control  will  enter  a modified  version  of 
portions  of  a program  by  Dennis  Phillips  and  Eric  Smith,  where  these  landmarks 
and  satellite  orbit  parameter  data  will  be  used  in  calculating  satellite 
attitude  parameters  necessary  for  the  final  transformations . Final  program 
output  will  include  time  sequenced  transformed  data  blocks  and  predicted 
parameter  values.  Expansion  of  output  capabilities  to  rough  prediction  of 
future  satellite  parameters  will  be  considered  at  a later  time. 


23 


CURRENT  INTERESTS 


The  automatic  landmark  recognition  program  uses  a previously  sectored 
version  of  the  original  data  tape.  A user  defined  area  from  the  tape  is 
stored,  and  coastal  outline,  cloud  covered  areas,  and  specific  landmark 
locations  arc  calculated  using  temperature  differentials  as  decision  criteria. 

Those  criteria  are  not  extremely  stable;  time  of  day  and  seasonal  variation, 

, 

as  well  as  weather  conditions  greatly  affect  the  reliability  of  results  on 
any  specific  trial. 

To  ensure  maximal  accuracy,  observed*  prevailing  conditions  and  anti- 
cipated approximate  coastal  outlines  of  each  landmark  site  are  considered 
in  developing  specific  recognition  criteria  for  each  similar  group  of  areas, 
e.,g.  different  approaches  are  applied  for  recognizing  an  island  than  for 
finding  a protrusion  or  indentation  in  an  approximately  vertical  or  hori- 
zontal coast.  Construction  of  a large  table  of  potential  landmarks  is 

important  to  guarantee  enough  usable  points  for  accurate  registration. 

IfJouM  6c 

Each  site  identified  by  approximate  location,  and  by  area  type  (indicating 
which  recognition  routine  to  use)/  Presently  work  is  being  done  to  analyze 
various  sets  of  landmark  data  in  both  the  visible  and  the  IR  to  eventually 
develop  a set  of  scan  routing  suited  for  most  landmarks  encountered.  A 
few  representative  data  cases  in  terms  of  temperature  extremes  are  being 
used  in  this  development.  However,  to  effectively  reduce  the  diurnal  and 
seasonal  uncertainty,  an  on-going  yearly  study  of  land-ocean  contrasts  at 
least  twice  daily  is  necessary. 

Prodipted  coastal  outlines  aYid  landmark  locations  are  visually  displayed 
on  a hard-copy  printer  where  manual  confirmation  of  acceptability  can  be 


m.ido  before  further  processing.  This  information  can  also  be  used  later  as 
a check  on  the  applicability  of  the  recognition  routine  being  used  at  each 
site,  and  updating  of  the  master  list  can  be  made  where  necessary.'  The 
capability  for  users  to  update  the  master  list  during  processing  will  be 
incorporated  into  the  routine  if  this  is  found  to  be  necessary  to  ensure 
the  acceptance  of  a large  enough  usable  data  set. 

«**  , 

Once  such  a set  is  identified,  the  approximate  location  tags  of  each 

site  allow  the  program  to  choose  those  landmarks  from  the  identified  set 

* 

that  minimize  later  error  in  the  calculation  of  transformation  parameters. 

For  this  decision,  landmark  separation,  to  be  maximized  for  points  equally 
distant  from  the  subsatellite  point,  and  distance  from  the  edge  of  the 
observed  earth  disk,  to  be  maximized  for  sets  of  points  with  equal  separation, 
are  considered-.  A final  set  of  landmarks  is  then  output  for  user  acceptance 
or  rejection. 

Tests  of  existing  software  have  shown  that  islands  and  previously 
defined  coastal  outlines  can  be  recognized  for  data 'collected  under  good 
conditions:  minimal  cloud  cover,  and  distinct . land /water  boundary  temperatures 
For  less-than-perfect  conditions,  further  testing  is  necessary  to  emphirically 
determine  optimal  coastal  recognition  criteria  adapted  to  anomalies  in 
regions  surrounding  prospective  landmark  sites. 

An  initial  attempt  was  made  to  obtain  necessary  transformation  para- 
meters from  a purely  geometric  standpoint.  The  transformation  from  earth 
reference  frame  to  photograph  line  and  element  was  qualitatively  correct, 
hut  accuracy  was  limited.  This  procedure  might  be  utilized  to  find  initial 
approximate  search  areas  to  be  processed  by  the  landmark  recognition  program 
described  above.  The  complexity  and  nature  of  the  inverse  transformation 


> I ' ‘ ' 

I , 

made  it  unstable  during  computer  testing,  and  the  occurance  of  negative 
radicals  caused  propram  termination  in  all  cases  attempted.  This  trans- 
formation was,  therefore,  insufficient  for  implementation. 

A more  sophisticated  approach  to  this  problem  is  a result  of  work 

done  by  Dennis  Phillips  and  Eric  Smith.  Portions  of  their  routine  are 

■ 

used  by  numerous  other  research  groups  in  this  field  with  adequate  results. 

r s 

A copy  of  software  in  NASA's  version  of  this  scheme  adapted  to  a mini- 
computer has  been  requested.  Two  main  procedures  from  this  program  are  to 
be  used  in  our  system.  The  first  utilizes  various  observed  orbit  parameters 
and  landmark  registration  results  to  compute  satellite  attitude  parameters. 

The  other  makes  use  of  these  values  in  the  block  data  transformation.  Final 

output  includes  user  requested  transformed  data  and  calculated  parameters. 

' • 

A long  range  satellite  position  predictive  technique  being  worked  on  by 
C.E.  Velez  and  his  group  at  NASA  Goddard  Space  Flight  Center  and  associates 

at  NOAA  NESS  has  not  had  as  good  success  but  their  progress  will  be  noted. 

\ 

s. 

I 
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FUTURE  PLANS 


Those  portions  of  Fhillip's  and  Smith's  program  applicable  to  our 

i 

proposed  system  will  be  modified  where  necessary  and  integrated  into  our 
program.  Studies  of  the  specific  computational  procedures  used  by  Smith 
will  be  used  to  determine  the  nature  of  an  optimal  decision  function  for 
the  acceptance  of  a final  landmark  set  from  those  points  successfully 
identified.  Tests  Involving  these  areas  will  be  completed  when  a sufficient 
table  of  identifiable  prospective  landmarks  has  been  compiled.  Also,  the 
possibility  of  incorporating  a cross -correlation  method  of  scanning,  the 
software  for  Which  was  provided  by  M.  Crowe  of  NOAA  NESS,  to  reduce  scanning 
error  will  be  investigated.  If  the  registration  scheme  developed  is  of 
reasonable  accuracy,  it  will  be  feasible  to  determine  cloud  top  height  and 
monitor  storm  development  in  a severe  storm  case  study. 
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APPENDIX  C 


Mathematical  Overview  of  D.  Phillips-E.  Smith  Report 


Toran  Hostbjor 


Sandra  K.  Weaver 


Rufus  E.  Bruce 
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m: 


[This  mathematical  overview  is  to  be  used  in  conjunction  with  the  D.  Phillips- 
E.  Smith  report  "Geosynchronous  Satellite  Navigation  Model"  referenced  in 
the  Bibliography  to  help  clarify  the  report.  All  terms  are  identified  in 
the  Phillips-Smith  paper,  as  are  the  equations  referenced.] 


Equation  (1)  gives  the  position  vector  of  the  satellite  in  the  earth 
coordinate  system. 

The  initial  vector  / N (t)  cos  ( \rr  (t~  tre^c  ) / P^)  \ 

l HCO  sin  J 

is  a vector  of  magnitude  H(t)  (the  satellite's  altitude)  lying  in  the 
earth's  equatorial  plane.  This  vector  is  operated  on  by  the  matrix 

A o O \ 

( o cos  Cx)  Sm(T)  J 

\0  -stnCz}  cos  Cl y 
which  rotates  the  vector  into  the  satellite's  orbital  plane. 

The  resultant  vector  is  then  operated  on  by  the  matrix 

/cos  (zrr  (t-  te^)/ Pc  ) s lnCs.Tr  o\ 

(-  sin  C 3/r  Ct-  t'fc')/ P,')  cos  farr  o 


which  rotates  the  vector  about  the  earth's  axis  by  an  angle  of  3 jj  & t 

~~K~~ 

Then  the  resultant  vector  is  operated  on  by  the  matrix 

(CosCxq  O -sinftiQc)  o\ 

sin  (pQc)  cos  feQc'')  O J 

O O 2_ 


which  rotates  the  vector  about  the  earth's  axis  through  the  longitudinal 
angle. 
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Equation  (3)  gives  the  position  vector  of  a landmark  measurement 
in  polar  spherical  coordinates 


K;  = 


fv  ces  &■  c.eS/\L 

rt" 

ry  sin  ey 


Where 

rC  = Radius  of  earth  at  landmark  i 
6>t-  = Latitude  of  landmark  i 
ft;  - Longitude  of  landmark  i 

— v 

Since  S(t)  is  the  position  vector 
of  the  satellite,  the  vector  C* 
to  the  landmark  i from  the  satellite 
may  be  written  as 


II;  U 


C-  = K-  - Sft) 


The  unit  vector  in  that  direction  is 
(Equation  6) 


K ~ 
i«,  - 


If  the  position  S(t)  of  the  satellite  is  known  for  a time  to,  then  the 
unit  vector  may  be  found  for  a later  time  t by  rotating  the  unit  vector 


at  t - tQ. 


K;  - 5VQ 
kv  - 5fol 


through  the  necessary 
angle. 


This  is  done  by  operating  on  the  vector  with  the  matrix 
/cos Ct  - 1„)  arr/fl,)  ~ ” (a  rr  (t / Pf) 

\ Snn  (a/r  ( t-  co%  (zrr  ft  -t„)/Ps  ) 


with  equation  (6)  as  the  result. 
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The  transformation  from  satellite  imagery  coordinates  (L,  E)  to  an 
earth  reference  frame  is  achieved  by  transforming  from  imagery  coordinates 
to  the  nominal  satellite  coordinate  system  and  then  transform  from  the 
nominal  satellite  coordinate  system  to  the  earth's  reference  system. 

Equation  (27)  gives  the  third  column  of  the  rotational  matrix  which 
performs  the  transformation  from  the  nominal  satellite  coordinate  system 
to  the  earth's  reference  frame.  The  resultant  vector  is  the  pointing  vector 
of  the  satellite  spin  axis  in  the  earth's  reference  frame.  "Tie  first 
column  of  the  rotational  matrix  is  obtained  by  the  Gram-Schmidt 
Orthogonalization  method.  We  have  the  vector  AROT3(t)  (column  3 of  the 
matrix)  and  we  have  the  satellite  pointing  vector  S(t).  Then  using  the 
Gram-Schmidt  method  we  can  get  the  first  column  of  the  matrix  AROT^ (t) 
and  equation  (30)  is  the  result: 


(0  = 


sc*y 

mu 

sti') 


. Azor  (i))  axotTCi) 

m\  J 


j 


This  is  equation  30.  The  second  column  of  the  matrix  is  a vector  perpendicular 
to  the  other  two  so  it  may  be  obtained  by  taking  the  cross  product  of 
AROTjft)  and  AROTj(t)  resulting  in  equation  (31) 

AZor^  (O  « ARor,  X *Rorx 

The  rotational  matrix  formed  by  these  three  vectors  is  given  in  equation  (32) . 
Equation  (33)  is  a rotational  matrix  which  corrects  for  the  misalignment 
between  the  camera  axis  and  the  satellite  spin  axis. 
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This  matrix  is  used  in  equation  (34)  to  get  the  pointing  vector  in  nominal 
satellite  coordinates.  This  vector  is  rotated  into  a pointing  vector  in 
earth  coordinates  using  the  inverse  of  the  matrix  in  equation  (32) . 

The  intersection  of  this  vector  with  the  earth's  surface  determines  an 
earth  coordinate  vector.  The  remainder  of  the  Phillips-Smith  report 
is  mathematically  straightforward. 
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MAIN 


DATF  = 76261 


16/02/4/ 


V G LJ!VFL  21, 

C MA  IN  PROGR  AH 

co^mcn/xl  ant/mi.  and » l day » icoof i ? 2 ) » pt  i m e 1 32 ) , xl ini  32),xflfi  *2) , 

1 XL  ATI  32)  , XLPNI  32 ) ,DlIN( 32) tOFLE (32) ,T IMELI 32  ) 

COMMON/ SC  ANR/I  ST  AN  , NUMSFN  , NOPC IN , TQTL  I N ,NFOL  I N , R ATI.  I N 

I , P ICLI N, TOTFLF, OcGFl F, P ADEL E,P I CFLF,EF, PITCH  , YAW  ,0.01 1 , $K  F W , w OTM I I 
2, p NT  MI 3, POTH21,ROTM2  3,ROTM31 ,ROTM33  ,RF ACT ,PF AS  I N ,S D* CO • POI R , PR AT 
COMMON/ GOAT A/PI, RCPCG,RF, A, B, A B, A $0,8  SO, ATMHGT , GRACON . EMFG A , SDL  5 1 D 
I » SH  A 1 1 YR » I HR 

CnMMON/NAVSlN/INAVtNAVNfLANDN»NIT,KIT#NnftB,KOAY#FL»FPf  ETtSPINO) 

1 ,P ASCEN,DFCLIN,SPINRA,TMPSCL,GTIM| 16), BET  A I 16  It  BOOT  I 16),NGAM( 16) 
COMMON/SYSCCM/  I T K , N L 
C 

PRINT  9 

9 FORMAT! //30X,22H  *+*  NAVIGATION  ***  /) 

ITK  = 1 

5 PRINT  810 

BIO  FORMAT! 22X, NAVIGATION  SOLUTION  ** ' ) 

IWO=  1 

13  GO  TO  1 14, 100,200) , TWO 

14  CALL  LOADMK 

25  CALL  LSORT(NLANC,PTIME,XLIN,XFLF,XLAT,XLONfI CODE ) 

16  NDA Y=L9AY 

MIO  = (NLANO+l  )/2 
HOUR  *PT I ME ( P I D ) 

CALL  GRT OR PI LO AY, HOUR  ) 

IB  SPINRA=100. 

TMP SC L=SP I AP A/3600000.0 
CALL  PRFPOS 
KB  ANO= 1 

CALL  SET  SCN I KB AND ) 

I JK  = 3 _ _ __ 

PRINT  999,IJK,PCLN,PICLIN 
I JK  = 9 

. PRINT  999, IJK,PCLN,PICLIN 

999  FORMAT! I 10,2E20.8) 

NA VN= 1 

__  PRINT  35,  NAVN 

35  FORMAT 1 1 8 X , * **  SPIN  ATTITUCE  SCLUTICN  NO.  • , 14, • **»,/2X, 

1 'DECLINATION  RT.  ASCEN.  CENTERLINE  LANDMARKS  SfARCH  ITER 

2 AT  I ONS ' I 

CALL  SPINSAXILANON  ,DEC  ,RAS  ,PCLN,MT  ,MIT  ) 

PRINT  40,  DEC , RAS « PCLN, LANON ,N I T ,MI T 
40  FORMAT (2F12.5,F6.0, • FIXED', 17,'  USEO',16,'  TOTAL ' , I 4, 'TURNS  • , 4X ) 
DFCL IN*OEC 
PASCFN*RAS 
PICLINaPCLN 
CALL  PRESAT 
CALL  RES l DU 

70  CALL  GAMCAL  _ 

GO  TO  707 
100  CONTINUE 

CALL  PREPOS 
CALL  PRESAT 
GO  TO  7B7 
200  CONTINUE 

GO  TO  787  * 

7B7  CONTINUE  1 

STOP 
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I 


n rv 


; LFVEL  21  FLTTME  OATE  = 76261 

$WROni)TINr:  FLTIKF(  1NT,  I DAY  * HniJR  ) 
niMc^^jpN  I N T ( 4 ) , MO AY (121 

DATA  MDAY/O, 31,59,90,  120,  151,  lfll*2l2*243*273»304»334/ 

tYR=TNT(l)/100 

MONaMnni INTI  1) , 100) 

IF ( MQN • L C • 1 2 ) GOTO  6 
JYP=(MON-l ) / 1 2 
I YO=I YR+JYR 
M0N=MfN-12* JYR 

6 IDA Y - 1 NT | 2 ) /l 00 

I DA  Y = ID  A Y *MD  A Y ( NON ) 

IF  1MDM.LT.3)  GO  TO  7 
IFIMPDI  1 YR,4) .EO.O)  1 DAY*  ID AY* l 

7 IDAY=I0AY»1000»U  YP-70) 

HOUR  = MODI  I NT | 2) , 100 ) 

FM I N=  I NT ( 3) /100 

SFC  = M0DIINTI3), 100) 

FSFC  = I NT  ( <► ) 

HOUR* HOUR* I FMIN*ISEC*FSEC/10000«0)/60,0)/60*0 
PRINT  20 

20  FORMAT! //2X  , ’ FLT IMF  : I DAY  , F M N, FSEC ,S EC , HOUR » ) 

PRINT  21,  1DAY,FMIN,FSEC, SEC, HOUR 

21  FORMAT  I /2X,I6,4F12»6) 

RETURN 

END 
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16/02/47 


k 


T IMOIF 


OA  TP  » 76261 


CVFL  21 

FUNCTION  T 1 MHl F ( IYROA1 .HOUR 1.IYROA2  » HOUR2 ) 


16/02/47 


IYl -vnn( IYRRAI /1000, 1001 
10  1=M0I)(  1 YR  P A 1 ( 1000) 

IFAP  l=(  IYl-1  >/4M 

01  = 365  < I Yl-l ) MFACUTCl-i 

I Y2=M0P(  1 YRRA2/1000,  l 00) 

1 1)2  = MOO  IIYRC  A2»  1000) 

I P AC  2= ( 1 Y2-1 ) /4  + 1 
02=165 M IY2-DH  FAC2M0  2-1 
Tl =1440.0*C1*60.0*H0URI 
T2= 1440.  O'*  0 2*60  • 0 ♦HOUR 2 
T 1 Mf)  TFs  T 2™T  l 
PRINT  10 

10  FORMAT!/ /2X»*TIM0IFs  I YRO A1 .HOUR  1 , 1 YR0A2 »H0UR2 tDl f 02 »T 1 » T2 , 
1 T I M C I F • I 

PR  I NT  11 »J YROAl f HOUR l « I YRO A2 ,HCUR2 » 0 1 , 02 * T l , T2 * T IM C I F 

11  FOR  MAT ( /2X  * I 10 » F 1 0. 4 » I 10. 3E 20. 4 , /2X» 3E 20. 4 ) 

RETURN 

END  " 


ao^ 


r.  L FVFl  ?\ 


I ROUND 


04  T£  ■ 76261 


16/02/47 


rtiNrnnN  ircuncu) 
f r f X I If  2*  3 

1 TRnufjn=x-0.5 
RGTUPN 

2 IPOUNO=0 
PCTIJPKJ 


3 IROUND*X«-0.5 
RETURN 
FNO 


on  — o noo  ooo  r*  o o 


IV  G LEVEL  21 


SATPOS 


DA  Tp  « 76261 


16/02/67 


SUB ROUT  IN':  S ATP0$  IN  AVOA Y , T IMP  » X , Y,  1 1 

CnMVPN/r.OATA/PI , R0PCG»  op  ,A,e,AR,  ASOf  rsq.  atmhgt,gracon,  emfga.solsio 

♦ tSHA.IYPtlHP 

Cn^MpN/SATOPP/  I0RR,  1X0 » XFR ♦ X I ! 3 ) t SI  AT , SLON,  SHOT , AR  I®  S . 

♦ t EOAY,FPHR, SEMINA, OCCC ON »GR 01 NCt FMANOM, PERHEL , ASNODE. 

♦ XMMC.SR0HF2,PX,PY,P/,QXtQY,07 

COMPUTE  MEAN  ANOMALY 

OIETIMsTI MOIF( IFCAY»EPHR,NAVCAYtTIME  » 

XMAN0M=XMMC»01FTIM 
ECANMlaXMANOM 
EPSILN-l.OF-8 

STIVE  FOP  ECCENTRIC 

00  2 I*  I » 20 
FCANOM=X  MANCM+CFCCFN* 

PRINT  20 

20  FORMAT I//2X, «SATPOS:XM 

PRINT  21 tXMMCtXMANOMt  E CANON 
21  FOP.MATI/2X*3E20.6) 

IF  I ARS(ECANOM-ECANMl).LT.EPSILN)  GO  TO  3 
ECANMl-ECANCM 

2 CONTINUE  . 

COMPUTE  CARTESIAN  COMPONETS 

3 XOMEGA-COS ( FCANOM  I-OECCEN 
Y0MFGA-SP0ME2*  S INI ECANOM) 

xs-xomega*px»yomega*qx  _ 

YS»XOMEGA*PY»YCMEGA*CY 
ZS«XOMEGA*PZ*YOMEGA*QZ 

PRJNT  30 

30  FORMAT l //2X, • SATPOS s XOPFGA, YCMFGA, XS t VS. ZS» I 

PRINT  31  t XCPEGA*  YCMEG  A*  XS*  YS » ZS _ 

31  FORMAT! /2X»5F15. 71 

ROTATE  TO  GFOGRAPMIC  COORDINATES 
a 

THR  * I HR 

DIFT1M«TIMOIFUVR,THR,NAVOAY,TIMEI 
RA«niFTIM*SriSID»PI/720.0C0*SHA 
RA$sAMOO! RA»2«0*Pl ) 

CR A-COS !R  AS ) . ' » 

SR  A*  SI  N!  R AS) 

PRINT  60 

60  F0RMAT(//2X,«SATPCS}DIFTIM,SHA,RAfRASvCRAtSRA«  ) 

PRINT  61,0!FT!M,SHA,RA»RAS'CRAfSRA 

61  FORMAT (/2X.6F20. 6) 

X*CRA*XS*  SRA*YS 
Y*«SRA*XS*CRA*YS 
7-7S 

XI  (1)»XS 
XI  ! 2 ) * YS 
X I ( 3 )■/ S 

SLAT -AT AN! 1/ SORT (X**2+V**2 ) J/RCPCG 
SL0N-ATAN2!  VtX)  /ROPOG  tr 


ANOMALY  1 V 

SIN! ECANM1 ) 

XMANOMtEC ANOM* ) 


4?  V. 


-VFL  PI  SATPQS  DA  TP  ■ 762M  16/02  /'♦  7 

• 

SHOT ■*  CQR  T ( X**2 ♦ Y* fc2*Z*'*'2l 
PRINT  10 

10  FOPMAT (2Xf’SATPCS:X»YfZ»Xl  (J)  WHERE  J*1 ,3, SLAT , SION . SHGT * ) 

PRINT  iltX«Y«Z.( XI( J) ,J«l ,31 .SLAT, SLON, SHGT 

11  FORMAT  </2X»<JF13«5) 

RETURN 

END 


PVPL  21  HALO  DATF  « 76261  16/02/47 

FUNCTION  II  ALO  ( X CFG  1 
INTFr,«=R*4  HALO,  U>EG, MIN, ISEC 
ZN^xncr, 

H I XIIFO.LT.O.O)  ZA*-ZN 
2 mnwN 
zn=  i nrr, 

7N  = t 7N-ZOH60.0 

PPINT  10, XOFG, IOEG,ZD,ZN 

10  FORMAT  I//2X,  • I L A L 0 : XDEC,  IO«=G,  ZO , ZN  * ,F  13. 6,  I 8,  2F  13. 61 
M1N=ZN 
ZMsMIN 

ZN* ( ZN*  ZM I ♦60.0+0.5 
I $5  C * 7 N _ 

HALO* 10000* IDF0M00*MIN*ISFC 
IFIXDEG.LT.O.O)  ILALC*-ILALO 
PRINT  20,MIN,ZM,ZN,I0EG,ILALC 


20  FORMAT I//2X, * IL ALO:  M INfZM, ZN, IOFG , ILALO *, I 8, 2F 13. 6, 21 8) 
RETURN 
ENO 


l 


LOAOMK 


OATF.  * 76261 


16/02  /4 7 


VEL  21 

SUBROUTINE  LOAOMK 
COMMON/ XL ANO/NL AND, LOAT, I CODE! 321 ,PTIME(32) ♦ XLI NI 32 1 , XELE C 32 » » 
l XL  AT  ( 32  I » X l C M 32 ) ,DL IN < 32  I , OEL E < 32  I , T IMEL < 32  I 
COMMON/BUFFFR/LPKNO,LMKIO,  IT  I PEI  6 I , LCOOE 

PENSION  JO  /Y  f 25 1 
INTEGER  v4  LAT , LCNG 
NLA  NO=LPKNC 
MO AY  = 32  700 
N*0 

LMK  IO  = -777 

IF  ILMKin.EG.-777)  GC  TO  200 
nr>  100  L — 1 » LPKNO 

READ  50,  LMKNn,LMKID,LAT,LCNG,XMP,XML,UTIME(J),J=l,4)  .LCOOE 

0 FORMAT(2I3,2I8,2F10.2,5I,5) 

PRINT  77,  LMKNO, LMKID,LAT, LONG, XMP, XML, (ITIPE(J),J*1, 4), LCOOE 

7 FORMAT!'  LPK*'  , 2 1 3, 2 18,  2F  10.2,  5 1 5) 

IFILCODE.LT. 01  GO  TO  100 
0 N-N41 

XLAT(N)=FLAIC(LAT)  _ _ 

XLONIN) *FLALO! LONG) 

XELE I N ) *XMP 
XLININ) »XPL 

CALL  FLTIME(ITIME,JOAY(N),PTIME(NI ) 

IF(JOAYIN).LT.  MOAY)  PDAY»  JO'AYIN  ) 

ICOOF IN) »LCC0E*100+LMK 1 0 „ _ 

00  CONTINUE  l 1 i 

NLAND*N 
LOATsMOAY 

DO  150  J*»l,NLANO  ' ' 

PT I me IJ) >PTIPE( J) 424.0*1 JOAY IJI-MCAY) 

50  CONTINUE _ i 

GO  TO  250  , 

00  LPKNO?  ' 

LMKIO-5 

1.C00E-0 

DO  5 1*1,4  . 

5 I T IMF  1 1 )*1  __ 

CALL  TESTMK 
PRINT  40  • 

40  FORMAT! //2X, 'LOAOMK :LMKNO,LMKI  0, ,LCODE  , I TIME  I II  WHERE 
PRINT  41 ,LMKNO ,UMKID»LC0CE»(  IT  IME (I) , I >1 , 4) 

41  FORMAT! /2X, 7181 

50  RETURN  

END 


G LFVFL  21 


TFSTMK 


DATE  * 76261 


16/02/A  7 


SUBROUTINE  TFSTMK 

COMMON/ XL  AND/NL  AND,  L CAY,  ICODE ( 3 2 1 ,PT I ME ( 32 1 , XL  I N(  32 1 , XCLF ( 3? ) , 
l XL  AT  (321  ,XLCM  12)  ,01  IN  (32)  ,DELF  (32  1 ♦ T IMFl  ( 32  I 
LPA Y»  A? 1 2 
Ml.  AND-* 9 

TLATMA. 65833333 
TL0N=-17.AA166667  , 

TIMP»12.0 
DO  20  I = 1 1 NIAND 

icpoe ( i >* l 

XL AT ( I ) *TLAT 
XLON( I ) =TLCN 

PTIMFU  )*TIME - 

TIMEsTf ME+0.5 
20  CONTINUE 

XL!N(l)*51A0 
XL1N(2)«50R8  I 

XL IN( 3 1 *50A  1 I I . ' 

XL  I N ( A ) *5002  j . . 

XL IN( 5) = A972 
XL  I N (61*4950 

XL1M  7)  = 4936  .._ . . ...  ' j 

XL  IN(  81=4930 

XL  IN (91=4934  l 

XELF  ( 1 1 *1 1A6^ ..... 

XELE ( 2 I * 1 IAAO  1 

XF  LF ( 3 ) *1 1A30 

XFLF( A) -11A20 

XELE (51*1 1A09 

XELE (6) *1 1A00 

XELE ( 71=11351  

XEL 5(81=11383 
XELE ( 9) = 1 1375 
• . PR INT  30 

3 0 FORMAT (//2X  »* TFSTMK »L DAY, NLAND,TL AT, TLON» ( XL  AT ( I > » XLON (I) ♦ PT IME ( I 1 
l»XLIN(I)*XELF(I) * WHERE  I *1 »9* 1 

PRINT  3l,LCAY,NLAN0,TLAT,TLCK» (XLATf I ) » XLON ( I ) » PTI ME ( I ) * XL IN( I I » 

1 XE LEI  I)  » I *1  *9) 

31  FORMAT(/2X»2l8t2El0.6f 1CF8. 2 » /2X, 1 5F8.2 , /2X #15F8.2 » /2X » 5F8 .2  1 
RETURN 
ENO  c 


* 


I 


-* 


43 


LFVFl  21 


LSORT 


OATF  * 76261 


16/02/67 


SlIRROUT  INF  LShRTINLtTJ, A,B,C»D.  !C  ) 

nivrNsinN  Tj(iitMiife(i),C(ii.ciiMcm 

L=  l 

10  TK=TJ?L) 

MV=L 

no  30  K=ltM. 

IFITJU1.GE.TK)  GO  TO  30 
TK*T J( K) 

MV  = K 

30  CONTINUE 
TS=TJ(U 
5A-MU 

SRsfMU  

SC*C(L) 
sn*n(Ll 
IS-  IC  IL ) 
tjiu.tk 

AIL  )*A  I MV  ) 

MLI-BIMV1  

C I L ) *CI MVl 
OIL  1*01  MV) 

ICIU-ICIKVI  . . ....  . 

TJfMVI-TS 

A(MV1*S A 

BIMV1-S8  

ClMV)»SC 

0|HV)»SD 

ICIMVI-IS 

L«L«1 

IFIL.LT. Nil  GO  TC  10 
77  RETURN  _ _ _ 

END 


r.  iFvri  ?\ 


PPFPOS 


SUBROUTINE  PRFPCS 

COMMON/ Gl)  AT  A /PI  ,PnPDG,R ,A,fl,AB,ASQ,BSQ, ATMMGT,GPACON,FMPGA,SOLS  10 
i ,sha, iyr, h-r 

COr'Mf'N/SATPRn/inRB.IXO.XHP.XSO)  , SL AT, SLON, SHOT, ARIFS, 

1 !FDAY,«;PHO,SFMIMA,OECCFN,ORBINC,XMANOM,PERHEL,ASNOOF, 

2 XVMC,SROVE?,PX,PY,PZ,CX,CY,OZ 
PI s3. 1A15S245 
F0P0G=Pt/lB0.0 

R* 6 3 71.221 
A=6378. 388 
R=6356.912 
AR*4^B 

ASQ= A*+  2 „ .... 

BS0=B**2 

ATMHGT=0.0 

SOLS  10= 1.00273791 

FMFGA=PI^SCISID/12.0 

SHA*  100. 26A67 

SH A» PDPDG^SP A . ...  . 

IYR=AOOl  .. 

IHR=0  1 ' 

GR  AC0N=0. 07A3657A 
RE*A  ' ■ 

XMMC=GRACCN*RP*  SGRTIRE/SEMIMAI/SEMIMA 
SB0MF2*SQRT  1 1 . O-OECCEN I + SCRT (1.0*0 ECCEN) 
n=PDP0G*0RBINC 
SOSINIO) 

CO=COSIO) 

P=POPDG*PERHEL 

SP=SIN(P)*$FMIMA 

CP=CO$(P)*SEMIMA  . Xr  * i 

S*RDP0G*ASN00F  , 

SA*S INI S ) 1 

CAacnsis) 

PX=CP*C A-SP*SA*C0  , 1 . 

PY  =CP*S A*SP*CA*CC 

PZ»SP*S0  ...  _ , 

CX*-SP*CA-CP*SA^CO 

QY=«$P*  SA*CP*CA*CC  ' 

QZ=CP*S0 
PRINT  10 

FORMAT! /2XV^PRBP0S:EMEGA vXPPCtSRCPEZf C*PfStPX«PYtPZiQX(QY«OZ* I 
PRINT  20,EMEGA,XHMCtSR0ME2fCtP,StPXtPY*PZ.QX»0Yt0Z  .... 

FORMAT I/2X ,3F12.8,AF20.At/2X»5E20.A> 

RETURN 

ENO  . _ . 


in 


PATE  » 76261 


16/02/67 


VFL  ?1 

FUNfTIPN  LO(IYR) 

I.P'C 

m MPP1  l YR  »4 ) .FQ.O)  ID*  1 

RF  Tl/RN 

END 


f 


G IFVFL  21 


PRFSAT 


DATE  ■ 76261 


16/02/67 


10 


11 


StlHPPUT  I NE  PRFSAT 

CfMMrN/SC/JNR/tSCAN  ,NUMS 

1 .PICLT N, THTRLFtOFOPLE ,R AOFLF , P I CF LF 
2,RnTMlJtRnTN2lfP0TM23,RnTM1lfRmy3^, 
errMPN/GOATA/P  I.RCPDGf  RF»  A,  6,  Af>,  A$Q, 
l * SM  A 1 1 VR  » IHR 

CnMPnN/NAVSLN/INAV,NAVN,LANON.NIT,MT 
1 »RA$CCN»0FCLIN»SPINRA»TMPSCL»GTIM(16 
PFC=OFCLIN*ROPOG 

sincfc*sinh:fci 

cnsDFC=ccs(rpc» 

i RAS=RA$CFN*ROPOG 

S I MPAS*  $ I N (R  AS  I 

COSRAS*COS(RASJ 

SP  INAX»C  OSOFC*COSR AS 

SPINAY»CCSOFC*SlKRAS 

SPINC  II -SPINAX 

SPINAZ=$INCEC 

SPINI2) aSPIAAY 

SPIN( 3)»SPTNAZ 

CP!TCH=RDPCG*PITCH 

CYAW=RDPnG*YAW 

CR  OL  L»ROPOG ♦ROLL 

PSKFW=ATAN2(SKEW,RA0L IN/PACELE 1 

PRINT  5 

FCFMATC//2X, 'PRESAT:  DFC  ,RAS»TMPSCL » 
PRINT  6fOEC«RAS»  TMPSCL »PS KEW ,S P I NAX» 
F0FMAT(/2Xt  7F17.EJ 
STP-SINICPITCH) 

CTP*CnS(CPITCHl 

STY*SIN(CYAW-PSKEW) 

CTY*COSICYAW-PSKEW» 

STR*SIN( CPOLL) 

CTRaf.OS  (CROLL) 

ROTMll*CTR*CTP 

R0TM13=STY*STR*CTP«-CTY*STP 

R0TM21--STR  _ 

' R0TM23*  STY*C  TR 
R0TM31=~>CTR*$TP 
R0TM33»CTY*CTP-STY*STR*STP 
RFACT-R0TM3  1**2*R0TM33**2 
ROASINa ATAN2CROTH3l,ROTM33 I 
PRINT  10 

FORMAT C 2X t * PRE SAT t SPIN ( 1 I *HEPE  I«1 
lROTM2l,ROTR23,RCTR31,ROTM33» I 
PR  INTI  l»  ISPINUI  »I«l,3»,R0TPiltR0TRl 
IROTM33 

FORMAT (2Xf9F13*6l 

RETURN  _ 

FNO 


PN.NOPCLN, TOTLIN.OFGL IN.RADL IN 
,FF, PITCH  ,YAW,POLL.SKPW,«>f)TMl  1 
RFACT ,RCASIN,SOt CDtPOI RtPPAT 
BSQ, ATMHGT,GRACON,FMFGAf  SOLSIO 

T,NORR,NnAY,ELfcPf ET,SPIN(3» 

It  BETA! 161 , BOOT! 16  I »NGAM( 161 


PSKE W » SPINAXfSPINAVtSPINAZ* I 
SPINAY t SPIN  AZ 


1 3 tROTMl lfR0TM13t 
3t  R0TM21 tROTM23t R0TM31 , 


1 | 
I 


47 


\ 


G LEVFL  21 


r,f:  THRU 


DATE  = 76261 


16/02/67 


SUBROUTINE  GETnpR(LDAYtHnUR  I 
r.r'VMHN/'MJI  F F» /KBAN0,MfMY 

cr^MrN/sATnwp/inRn, ixd,xhr,xs{ 3),slat,slon, shot, ar i e ^ • 

1 TFHAY.r phr, SEMIMA.OCCCEN.OREINC.XMANCM, perhel, asnode, 

2XMMC  ,SR’)MF2,PX,PY,PZ  ,QX,QY,QZ 
IfjRn  = -l 

1H10RBI  10»70»30 
10  !FDA Y=62 16 
FP  HP =0.0 
SEMIMA=62168.86 
CFCCFN=0. 001207 
HP  R I NC  =1 . 920 
XMANOM* 181.235 
PFPHEL=267.316 
AS NODE =158. 189 
GO  TO  70 

30  RF  AO  33, IEOAY,FPHR  , S FM IMA , OECC EM* OP RI NC , XM ANOM, PERHEL , A SNODE 
33  FORMAT! I 6*7F15.6) 

ifoay=mday 

PP I NT  65 

65  FORMAT! //2X,*GET0R8s  I EOAY ,FPHR , SEMI MA ,0ECCFN,0RBI NC , XMANOM, PERHEL 
1, ASNODE* I 

70  PR  T NT  80  * I FD AY , E PHP , $ EM  I M A ,C ECC EN, CR B INC , XMANOM,  PERHEL , ASNODE 

CALL  EPOCH! I EDAY , EPMR * SEM I MA ,OECCEN , XMANOM 1 

PR  I NT 80,  If  CAY, EP HR, SEM IMA ,OECCEN, ORB I NC , XMANOM, PERHEL ♦ ASNODE 

80  F0PMATI5X,* ORBIT* ,I6,7F15.6) 

RETURN 

END 


rvFL  ?l 


SFTSCN 


OATf  = 7<S?6l 


16/02/47 


SlJRROUTINF  SETSCNIKRANOI 

CCM^TM/SC  ANR/  ISC  AM  , NUMSFN  ,Nnpr  i_N , TDTL  1 N , DFGL  I M , R AIM.  I M 

l , flCLTM, TO  TEL  E ♦ OEGr IF  ,R AnfLE » P I CFLF , TF, P I T CH  , Y Aw .POLL , SK FW , »0 TM l 1 
2,enTMt  l,ROTR2l.ROTM23»ROTM31,ROTM33tRFACT,ROASIN,SD,CD,POIR  ,PRAT 
C OMMOM/qo AT  A/ P 1 1 ROPOG » R E » A * P» A R, ASG.BSQ, ATMHGT , GR ACON, EMFG A , SOL S I 0 

1 ,$ha,iyr,ihr 

N0PCLN=0 
PREP  AT  =0.0 

ppeniR=o.o 

PIT  CH=  0 • 0 
YAW=0.0 
onLL  =0.0 

SKFWrO.O  . 

PO I D = 0«  0 
PRAT =0.0 
sn=5TNtPniR) 

Cn=COSIPDIR)  , 

IF  IKPANC.EQ.2)  GO  TO  30 

NUWSFN  = fl  , __  , 

GO  TO  A 5 
NUMSFN=2 
GO  Tf)  45 
I SFN$nB=NUMSFN 

TOTlIN*182l.O*SFNSCP  , 

DEGL IN=  20.0 

TCTFLF*  lqi  l.O^SFNSOR 

DEGE IE  = 18.375 

PICL  IN=  ( TOTLIN+l.O) /2.0 

PAnilNaRDPOG^CFGLIN/ITQTLIN-I.O) 

RAOELE=RDP0G*nEGELE/<T0TELE-l.0) 

P ICELE=I L .0*T0TELP)/2.0 
EF=R  AOEIF/(2.0*F1) 

PRINT  53t  KBANO,  SENSOR  ,TOTL IN, DEGL  IN,TOTELE ,DEGELE, PICL IN, SO 
FORMAT ( / 16H  SCAN  CONSTANTS.  ,15  , 2F  13.  I , F 1 3.  4 , F 1 3.  I ,F  1 3. 4, 

1 2F13.1J 
I JK  a 1 

PRINT  999, I JK,PCLN,P I CL  IN 

FORMAT  ( • PCLAtPICUN  PRINT  NO.  • , 1 7, 2 ( 1 PG2Q  . 10  )) 

RETURN 

END 


dSL 


V G LEVEL  ?1 


epoch 


DA  TF 


76261 


16/02/67 


RlfO&nilTINF  FPncm  I FT  1MY.EPHR  , SEM ! MA  .GECCEN  , XMF ANA  1 

PI "3. 1 61 592C5 

ROPDGsP 1/100.0 

Re -637H.  3«6 

GRACON- 0.0 1636576 

Xmmc=GP  ACTN*  SQP  T!RE/SFMIMA)«'t3 

xMANnM=wnpoo*xvF  ana 

TlME=(XMANrF-CECCFN*S INJXMANOM) l/(60.0*XMMC  ) 

PRINT  20 

20  Fr»MAT< //2X»  'EPOCH:  XMMC  , XRANQM , XME  ANA , T I ME  ' ) 

POINT  21  tXPRC,XRANCMf XMEANA.TIME 

21  FORMAT! /2X, 6FI 5. 71 
TIME1=FPHR 
TIMr=TIMFl-TIMF 

inA  Y = T I mp /26 . 0 

PRINT  30,TIPE, ICAY 

30  F n R M A T I /2X* 'EPOCH!  T IME ,IDAY* ,/2X,F12.6,l6) 

IF! T IMF.LT.0.0)  I C AY=  I DAY™  1 
TI ME=TI me-26.O*I0AY 
<t  FPHP=TI ME 
XMRANA=0.0 

PRINT  6 0 f E PHP , XMEANA 

60  FORMAT! /2X, 'EPOCH:  EPHR  , XMEANA * , 2X  ,2F 12. 6) 

IF ( I DAY  * FO.O ) GO  TO  12 
JYFAR=M0D! IFTIMY/1000,100) 

JDAY=MOD( IFTIMY, 10001 


7 

12 


JO  A Y = JDAY ♦ I C AY 

IF! JDAY.LT.  II  GO  TO  5 

JTGT=365+ID!JYFAR) 

IF! JDAY.GT. JTOTI  GC  TO  < 

GO  TO  7 

JYF AP=j YEAR"l 

JDAY=365*LD< JYEAR I ♦ JO AY 

GO  TO  7 

JYEApsJYFAR+l 

JOA Y* JDA Y" J70T 

I FT  IH  Y*  1000  ♦ JYEAR  4- JOAY 

RFTURN  " 

ENO 


I 


' s 


!..  i 


mu. 


so 


— ■ — : — -rr 


\ 


r,  irvFi  21 


SPINA* 


DATF  * 76261 


16/02/6  7 


30  FORMAT  ( / / 2 X » ' SPINAXs  XFAC  T , UX,  IJY.tJZ  , YL  IN  • , /2X , 5F  13.  51 
MNI IN=SIN(YIIN) 

CflSI.  IN=CDS( YLIN) 

nn  )*oii  i+ux**2 

0(2  ) =r>(,2  ) ♦LX-UY 
1 r(  1 )-n(  3 ) MJX  *IJ  7 
• t * l -i*l  ( * 1 1 V * * 2 
'.*(  M -n(  *,)  ii)t*  ij/ 

>'((. ) = I*  ( 6 ) « U / * * 2 

ni  n =o<  n hjx^ccslin 
CIO)  =r»(  8 I *UX  *SI  NLI  N 
n(9>=r>(9)fUY*ro.$L  in 

n(  1 0 ) =r><  10)  fUY+SINLIN  

n< 1 1 »«o( 1 1 )«iiz*coslin 

H(  12)  =0  ( 12J+U7  + SIM.IN  , ' i 

D( 1 3 ) =D I 13) ♦C0SIIN**2 
nil6)*0(l4)+SINl  IN**2 
n< 15)=D( 15)+C0SLIM$IMIN 
P° INT50 

50  PflBMAT(2X*,SP!KAX!Xl»X2»X3*C(J)  WHERE  J»l,15') 

PRINT  5l,Xi,X2,X3t(0( J) , J»i,15) 

51  FOPMAT (2X,9Fl2.flt/2X,9F12.8  ) 

3 CONTINUE 

IF  ( NUMSPN.EC. 0 ) RETURN 
IF  (MUMSPN.EQ.11G0TO13 
NITLIM*5000 
HITt  IM«5CiOO 
XL  AV =0.01 
ZFRQT*l.F-7 
Sl=0.0 
S2*PI/2.0 
S3-PI/2.0 
PARM1«0.5 
. PARP2-0.75 

PAPM3— 0.5 
N«1 
!•# 
j»o 

4 r» i*i 

5 IF! XLAM.LT.ZERCT1G0  TO  11 
SA*DS I N<  SI)  * 

CA=rcos (si ) 

$B»0$IN( S2 ) 

CB=PCO$<  S2  ) 

SP-0SINIS3) 

CP-OCOSI S3) 

GA  = C(l>*CA*SA*C(2  )*C8*I2  .0*CA**2-1.0)  + 0(  3)  * I 2.  O-CA*^  2-1 . 0)  ♦ SB-0  U) 
l*CA*SA/*CB»*2-O(5>  + 2.O*SA+CAfcSB<CR-D(6)*SA*CA*$R'»*’'2-0m*CA*CP*O(  B) 
2*CA*SP  +0  1 9 >*SA*CR*CP-D(  10  > * $ A*C B + SP  + OI  1 1 )* SA*SB*CP -0 1 12 1 *S A* SB*SP 
GB*-0(2  )+CA*SA*SB*0(3  )*S  A-*CA  ‘Ce-0(4  I *CA  +-*2 *CB *SB+0{  5 ) *C A++2+ C 2.  0*C 
IB*-  2-1.0) +01  6l*CA  + »2*SH+Cn«-l)t*))*CA«SB*CP">n(  10 ) ♦CA*SB*S  P-0 1 1 1 )+CA*C 
2B*CP+n(12)  + CA*CMSP 
IF (N0PCLN.EC.0.CP.NUMSPN.FQ.2)  GG  TC  6 

GP*0(7 )*SAR$P+DI 6 )TSA*CP+0I9)*CA*CB*SP+0I 10I*CARC8*C P+0 I l l ) «CA*S8* 
1 SP+P ( 12  )+C A+SB^CP-DI  13)*SP*CP+IM14)*SP+CP+D(  15 )« I 2.0*SP**2-1 . 0) 
PRINT  60 

60  F0RMAT(2X,*SP!NAX:GA,GB.GP*) 

PRINT  61  , GA.GB.GF 

52 


l rvri_  ?\ 


SPINAX 


OATE  r 76261 


16/02/47 


M * u«.H  (/X  .3F12.81 
Gil  TO  7 

6 rp=n.o 

7 r,M  =fiS0R  T (f.  A**2*CB**2  *GP**2  ) 
ga=ga/gm 

gb=gr/gm 
gp=gp/gm 
GO  T0( 8 » 10 1 ,N 
B N=2 
9 G1=GA 
C2  = GB 
G3=GP 
XA*S  1 
XR*S2 
XP  =S3 

S 1 = XA-Xl AMV  GA 

S2=XB-XL  AM^Ge 

$3=XP-XLAM*GP 

IFf  I .EQ.NITLIM)  GO  TOH 

GO  TO  4 

10  COSA=G  l*G  A+G2'*G8*G3*GP 

XIAM»XLAMMC0SA*PARM1*PARM2) 

IF<C0SA.GT.PARP3)  GC  TC  9 

S l*XA 

S2*  XB 

S3=XP 

J»J*1 

1F(J.F0.MITL IMIGO  TO  11 
GO  TO  5 
11  PRINT  70 

70  F0RPATI2X,*SPINAX:C0SA,XLAP» » 

PRINT  71,  COSA ,XIAM 

71  F0RKAT(2X,2F15.71 

• MIT  » I 

MIT-J 

SPAX1*  OS  IN ( S 1 ) 

SP AX2 *OCOT<  S1)*0CCS(S2)  _ 

S<»AX3»0C0SIS11*0$IN(S21 

0FC*90.0«ATAN(S0RT| SP AX  1*  *2  + SPAX2**2 1/SPAX3 1 /R  OPDG 
PAS*0. 0 

IF( SPAX3.GT. 6. 999999991  GO  TO  12 
RAS*ATAN2(SPAX2,SFAX1)/R0PCG 

12  PClN*PCLN-(S3-PI/2.0)/RACLIN 
PRINT  80 

00  F0PHATI2X,* SPINAX:  CEC ,P AS ,PCtN* 1 
PRINT  81,  DEC ,RAS,PCLN 
81  FORMAT (2X,3F15«71 

13  CONTINUE 
RETURN 
END 


:VEL  21 


satfar 


OATF 


76?#,  I 


16/02  /'♦  t 


SUFFOllTfNP  SATFARlPtCTlM,KLTN,XELE.,XLAT,XLON,  ! TYP  F,  NFRR  , BET  A IN,  T 
l DOT, ATFRAC I 

fFATF  Rr>G:i350,6)5ATFAP.FIN/NV 

SATEAP  COMPUTES  SATELLITE  rjiQW  * EARTH  COOR  * FARTH  EDGES  ♦ SUB  POINTS 
TIOI  IS  HEFINFC  TO  PE  GREENWICH  HOUR  0 OF  NAVIGATION 
IATITVJOF  RANGES  FROM  490TC-90  SCUTH 

LONGITUDE  RANGES  FROM  ♦10OTO-1BO  WEST 
*♦*•.*• 

INPUT  PARAMETERS 

PICITIM  = PICTURE  START  TIME  ( POUR  FROM  TIOI  I 

XL  I N * SATELLITE  COCROIINATE  I LINE  I ' ' 

XFLE  * SATELLITE  COORDINATE  I ELEMENT  I - 
XL A T = FARTH  COCRCINATE  I OEGREES  LATITUDE  ) 

XLON  » EARTH  COORDINATE  ) DEGREES  LONGITUDE  I 

ITYPE  = 1 FOR  SATELLITE  COORDINATE  TO  EARTH  COORDINATE  TRANSFROM 
= 2 FCR  EARTH  COORCINATF  TO  SATELLITE  COHO IN AT E TRANSFORM 
» 3 ECR  LEFT-RIGHT  OBLATE  EARTH  EOGE  C XL AT  = LEFT  , XLON  - RIGHT  » 
» 6 FCR  SUP-PCINT  ( XLIN-LINE  , XELE=ELF.  ,XLAT«LAT.  , XLON«LON.  I 
: * 5 FOR  ROTATION  ANGLE  I XlIN  * ROTATION  ANGLE  I 

NEPPs  ERROR  FLAG  (=0  FOR  NORMAL  PETURN,  « 2 THRU  9 FOB  ERRORS! 

RFTAIN  * BETA  ANGLE  AT  TIOI  I FL  FMFNTS  1 

PE TDOT  a RATE  CF  CHANGE  OF  BETA  I FLEMENTS  PFR  HOUR  I 

AT F R AC  * CLOUD  HEIGHT  COEFFICIENT  I RANGES  FROM  0 TO  l | 

; ITFR  e I TFPATION  COUNT 

GAMMA  = beta  ANGLE  AT  SAMPLE  POINT  TIME  I RADIANS  I 
SAMTIM  a SAMPLE  POINT  TIME  I HOURS  FROM  TIOI  ) 

C.T.  MOTTERSHEAO/CSC 
713  OCT  1975 


CCHHCN/GOATA/PI, RCPCC.R  , A, B, AB , ASQ, BSQ, ATMHGT , GRACON, FMFGA , SOL S I D 
I , SHA  » I YR  » I HR 

• COMMON/ SOANR/I  SCAN  ,NUMSEN  ,NOPC  IN , TCTL I fl.DFGL  T N ,R  AOL I N 

It  POLN  tTOTELEt  CEGELE,RADELE,PICELEt  FF,P !TCH,VAW,ROLLtSKEW,R0TMllt 
2POTM13, R0TM21,R0TM23,P0TM31 ,R0TM33 t REACT, RO AS  IN tS Pt CCtPO  IR t PRAT 
COMMON/ N4VSLN/ I NAV,NAVN,L ANON, NIT, MI T,NORB, NOA V, FL ♦ FP, ET, S P INI3I , 

I RASCEN,DFCLIN,SPINRA,TMPSCL,GTIMI 16  I, BET A I 16),B00T| 16 1, NG AMI  16) 
COMMON/ S A TVEC/X SAT ,YSAT,2SAT,  HEIGHT ,XVECl «XVEC2 ,XVEC3»THET  A, 
1AR0T11,  ARQT12, AR0T13 , AR0T2 1 , AROT22 , AROT23, AROT 31, AROT32 1 AR0T33 
LIMIT-3  * 

NERR-0 

BRANCH  CN  TRANSFORM  TYPE 
GO  TO  112,9,12,9,9),  ITYPE 

INITALT2E  FOR  LATLCNG  OR  SUBSATELL ITF  POINT  CALCULATION 

9 XLINaPC.LN 

XFL  F*  PI  CEL  E 
I TER -0 

COMPUTE  TIMF  OEPENOENTS  PARAMETERS 

12  SAMTlMaVTIMFIPICTlM,XL  IN  ,XFLF  ) 

GAMMA aRAOELE*! BET AIN* BE  TOOT "SAMTIM) 

CALL  SATC0RIN0AY,SAMT1M,GAMMA,SPINI 

54 


1 6/ 02 /A  7 


DATE  = 76261 


P UNT  in, SAMT|M,PJCTlwtXLINt XFLF,OAMMA,RFTAIN,BFTOOTf SPIN 
10  FOPMAT  I //2X»*  SATFAR  : S AM  T I " , P I CT  IM,  XI.  IN , XFL F , G AMMA , B FT A IN , BF TOOT f 
ISPlN't /2X,AE20.fc,/2X,AF20.6> 


CHFCK  FOR  TRANSFORM  OIPFCTIfN 


TRANSFORM  FROM  SATFLL I TF  COORDINATES  TO  EARTH  COORDINATES 
**  **  * 

C At  L LATICM  ATFRAC»XLIN,XFlEtXlAT,XLCN,NERR  I 

PR INT  TO  1 I - ' 

FORMAT (2X, »SE:  XLATfXLCNf NERP* > 

PRINT  3l,XLAT,XL0N,NFPR 
FOPMATI 2X,2F15.7, IBI 
GO  TO  77 


TRANSFORM  FROM  FAPTH  COORCINATES  TO  SATELLITE  COORDINATES 


TTFP=ITER*1 

I F < ITFR.GT.LIMITIGOTO  18 
TF<  ITFR  .GT.l  I GOTO  17 

t.**.  + * 

COMPUTE  EARTH  COORDINATE  VECTOR 

CALL  VCLALOI ATFP AC t XL AT  , X LON f XE » YF, Z5, NERR ) 
PRINT  AO 

FCPMAT (2Xt • EAPTH  COnR . VEC.J  XFtVE, ZE» NERR  • I 
PRINT  A 1 , XE  »YE  * ZE  » NERR 
FQRMATC2X,3E20.A,I8)  _ 

IFtNFRR.GT.0I  GOTO  77 
CALL  LINELr<XE»YEtZE»XLIN»XFLE) 

PRINT  50 

F0RMATI2X, • LlNELEl : XL1N*XELE') 

PRINT  5 1 1 XL  IN,  XELE 
FORMAT (2X,2F15.7I 
GO  TO  12 


CHECK  IF  POINT  IS  OFF  FRAME  ANO  IF  SO  SET  ERROR  FLAG 


I Ft  XLIN.LT  .1.0  .OR. XL I N. GT • TOTL  TNI  NERR-A 
IFtXELE.LT.l.O.CR.XELE.GT.TOTELEI  NERR-5 
GO  TO  77 


***** 

EARTH  EDGE  COMPUTATION 

COMPUTE  POINTING  VECTOR  IN  SATFLUTE  COORDINATE  SYSTFM  AT  FLFMENT  0 


SUBROUTINE  HORIZON  NOT  USED 
70  CONTINUE 
GO  TO  77 


SUB-SATELLITE  POINT  COPPUTATiCN 


uuo 


SATFAR 


DATE  ■ 76261 


16/02  /A  7 


CALL  LINFIF<0.0,0.0,0.0,XLIN,XFIE) 

P»  TNT  60 

60  FORMAT(2X »• LINEIE2:  XLIN.XFLE'I 
PRINT  61,  XL  I N » XELE 

61  FORMAT, (2X.2F15. 7) 

' C.O  TO  12 

COMPUTE  SUn-SATELLI TF  POINT  FROM  SATELLITE  POSITION  VECTOR 

23  XLAT=ATANCXVEC3/SCRT( XVEC1* >2*XVEC2T*2 ll/POPOG 
XL  FN=  AT  AN2  IXVFC2 , XV f C 1 ) /ROPOG 
PRINT  70 

70  FOPMATI 2X. • SUB-SAT  PT.:  XLAT,XLCN»)  

PRINT  71,  XL  AT , XLON 

71  FORMAT(2X,2F15.7»  • 1 

IF!  ITYPE  .EQ.6I  GO  TO  77 

77  RETURN 
END 


n no  r»  o o r»oo  r»oo  onor>or»c>r>r»noo 


SATCOR 


OATf  = 762M 


16/02/A7 


’ L EV  Ft  21 

SURPOHTINF  SAT  CC°  (NO  Ay,  S A”T  JM,  GAMMA,  SP  IN  ) 

SATOR  USCS  THF  NAVIGATION  SOIUTICN  ( OR  R I T , S P INATT  I TUOE  . GAMM  A ) 

TO  CAL  GUI  AT  F THF  SATELLITE  COORDINA  TF  ROTATION  MATRIX  AROT 
AT  THF  SARPLF  tjmf,  AND  STORE  IT  IN  THE  COMMON/S ATV EC/ 

SG  s SIN  ( GAMMA  » 
fC.  = COS  ( GAMMA) 

CARROT  = EARTH  ROTATION  FROM  TIO)  IPAOIANS  ) 

ST  « SIN  C EARPOT  ) 

CT  * COS  I EAPROT  ) 

C.T.  MOT TER SHE  AO/CSC 
10  OCT.  1975 

COMMON/GOATA/PI  fRDPr)GtR*A»5fA8*ASC*BSO»  ATMHGT,  GR  ACON»EMEGA ,SOLS I D 
* »SHA  , I YP* IHR 

COMMON/ SATVEC/XSATfYSAT»ZSAT*HFIGHTt  XVEC 1 , XVEC  2 » XVFC3, THET  A* 

• ARCTll,  AROT  12, AROT 13 , AROT 2 l, AR0T22 , A»OT23 , AROT 3 l, AROT32 , AROT33 
OIMENSICN  SFINI3) 

SG=SIN(GAMMA) 

CG»CCS I GAMMA) 

PRINT  30, GAMMA, SG,CG 

30  FOPMATI //2X, ' SATCOR : GAMMA , SG ,CG * , /2X, 3E20. 8 ) 

CARROT«PMFGA*SAMT IM 
iST*  S I N(  E'ARR OT) 

CT»COSIEAPROT) 

PFPOINT  SPIN  AXIS  AS  FUNCTION  CF  PRECESSION 

SP  A X 1 =SP I N ( 1 ) 

SPAX2=SPINI 2) 

S PAX3*S  P IN ( 3 ) 

CALL  PRFCESI SAMTIM,SPAXl,SPAX2,SPAX3) 

i 

• COMPUTE  DISPLACEMENT  VECTOf<  FROM  ORBIT  INFORMATION  C KEPLER  I AN  MOOEL  » 

CALL  SATPOSINDAY,SAMTIM,XSAT, YSAT.ZSAT) 

PRINT  50 

50  FORMAT! 2X,*PRECES  SATPOS:  SPAXl » SPAX2  »SPAX3  »X$AT  ,Y$AT»  ZSAT* ) 

PRINT  51,  SPAXl,SPAX2,SPAX3*XSATtYSAT,ZSAT 

51  FORMAT(2X,6F^5.7) 

COMPUTE  UNIT  POINTING  VECTOR 


HEICHT=SQRT(XSAT**2»YSAT**2AZSAT**2) 

XVECl-XSAT/HEIGHT 

XVEC2*Y SAT/HEIGHT 

XVEC 3* 2SAT /HEIGHT 

THETAsR /HEIGHT 


rOMPUTR  NOMINAL  SATELLITE  POSITION  ROTATIONAL  MATRIX 


AROT31»CT*SPAXl*ST*SPAX2 
AROT32^-ST*SPAXl*0T«  SPAX2 
AROT33»SPAX3 

COSA=XVFC  l*AROT3l*XVE02#AROT32*XVEC  3* AROT 3 3 
YVrn-XVrCl»COr,A*AROT3l 
YVFC2-XVFC2-COSA*ARCT32 
VVFC3»XVFC3-COSA*AROT33 


r»7 


t 

\ 


evrL 


$ ATC.PR 


0 ATF  = 76261 


16/02/47 


21 


YfiiOfr  -*t -OAcrT » vvFCl**  ?*  yvi  c?**  ?»vvrr. 3**2 | 

PfU.Mr  40 

40  POl’MAT  t ///X.'SATCCR  s FAR«CT,HriGHTtSAMT!M,cnSA,YVECl  fVVFC2,VVFC3f 
iVM'TR  • I 

PRINT  41  ,FAPRrT,HF!OHT,SAMT|M  ,COC A f YVEC l f YVEC2 t YVEC 3 , YNOR 

41  FORMAT! /2X.PF1J. 5)  ' 

A Pf T 1 1*  YVFC l vYNOR  1 , 

APOT12*  YVFC2t'YN0» 

APOTH  = YV*  C3*YN0R 

APnT21=AQC;T32*APCT13«APOT33*AROT12 
APnT22-ARnT33*-ARCTl  1 - AR CT 3 l * ARO T l 3 

ARPT?3t  AP0T31*  AROT 1 2- ARPT 32*ARO T l 1 1 

ARPT  sAROTl 1 1 | 

AROT  i l=APPT*CG“SG*ARRT2  1 
ARnT21-ARnT*SG*CG*AR0T2l 

ARC  T = AR0T1 2 

ARnT12=APnT«CG-SG*ARnT22 
ARnT22»ARPT  + SP.*r.G*APCT22 
AROT  = AP IT  1 3 

AROT 13=  AoQT*CG-SG*AROT22 
ARnT23=ARCT*SG«r.n*ARRT23 
PRINT  60 

60  FORMAT !2X» ' SAT COR  : AROT 1 1 , AROT 12, AROT 1 3, AR0T2 1 , ARQT22* AROT23 
l,ARnT31,AP0T32,AR0T33* ) 

PR  TNT  6l,AP0Tll, AROT 12 , AROT  13, ARCT2 1 » AROT22 t AROT23 1 AR0T3 1 » 

1 APOT32  t ARCT33 
61  FORMAT! 2X,9F13. 6) 

RETURN 

END 


LEVEL  21 


L4TL0N 


RATE  = 76261 


16/02/67 


C 

C 

C 


SUPROUT  I NE  LAHONI  ATFRAC,  XLlN,XFLrfXLATfXLnN,NFRR) 

CCMKrN/r.OATA/Pf  ,RrPnGtR,A.«t  AO,ASQ»HSQ,  ATMHr,T?r.7ACnN,EMFGA  f sniSIO 
* t SMA , l YR  » IHR 

CTM^ON/  >C  ANR/  1 SC  AN  , NUMSFN ,M(lPCl  N,  TOTl  I N *0FGL  \ N ,P  AOUN 

l .PICLIN.TnTFLE.nEnFLF.RAnFLF.PTCFLE.EFtPITCH.YAW.RnLL.SKFW.onTMU, 
♦ROTNM,  ROTP2  l,ROTM2  3,POTM3l  , ROTN33  ,RFACT  ,RCASl  N,S0,C0,  POl  » , PRAT 
rCMMrN/SATVrr/XSAT,YSATf2SAT,HFIGHTf XVCC1, XVFC 2 , X VEC 3, THF TA , 

■*  AR  nr  1 1 , ARTTl2t  AR0T13,AR0T2l»  ARCT22,  AR0T23,  AP0T31,  AROT32fAROT33 
YL I N*l XL  I N-PCLN ) *RAOL I N 
YFLC=(XCLE-PICELEI*RACELE 
SINLI N*  SI  Ml YLIN) 

COSL IN=COS  I YL I N I 
SINriFsSIN(YELE) 

COSFLE*CO$(YFLFI 

COMPUTE  POINTING  VECTOR  IN  SATELLITE  COORDINATE  SYSTEM  AT  ELEMENT  0 


C 

C 

C 


C 

C 

C 


C 

C 

C 


C 

C 

C 

c 

c 

c 


FLT=RCTM11*C0SL IN-RnTM13*S INLIN 
EMI»R0TM2l*C0SLIN-R0TP23*SIMIN 
ENI»ROTM3l*COSLlN-RnTM3  3««SINLIN 

ADJUST  POINTING  VECTOR  FOR  ELEMENT  COUNT 

FLI  sCOSELE^ELT+SINELF^ERI 
FMI»-SINELE6ELI»CCSELF*EMI 
FM»ENI 

COMPUTE  POINTING  VECTOR  IN  EARTH  COORDINATE  SYSTEM 

ELO=AROTll*FLI+ARCT2l*FNI*AROT3l*FNI 
FPC=AR0TI2*FLI-»AR0T22*EM  I4AF0T32*FNI 
FNn  = AR0T13*FLl*ARCT23>*FPI-»AR0T33*FNI 
. PRINT  10 

10  FORMAT ( //2X,'LATL0N:  YL  IN , YELE , EL  I ,EMI  , ENI ,FLI , FMI ,FNI ,ELO,FMO, 
1 CN0 • ) 

PRINT  1 1, YL  IN, YELE t EL  I, EM  I, ENI, FLI ,FMI ,FNI , FLO,EMO,ENO 

11  F0PMAT(/2X,11F10.3I  

ADJUST  FOR  CRLATENESS  OF  EARTH  SPHERE  AND  CLOUD  HEIGHT 

9 

cldhgt*atfrac*atmhgt 

AHGTSQa ( A*CLCHGT)**2 
BHGTSQ*! B*CICHGT>**2 

basq»rhgtsq/amgtsq 

CNFMSC-1 .O-EASQ 
AQ«PASQ6CN'PS0«ENC**2 

BQ»2.0*t I EL0*XSAT6EM0*YSAT )*baso+enq*zsat) 
CQ«<XSAT»*2*YSA'r**2)«pASC*7SAT**2-BHGTSQ 
RA0«BQ**2«6.0*Aq*CQ 

CHECK  IF  POIT  IS  CFF  EARTH  ANC  IF  SC  LFT  REJECTION  VALUES 
IFIRAD.LT. 1.0)  GO  TO  32 

FIND  POINT  ALONG  POINTING  VECTOR  INTERSECTING  FARTH  SURFACE 

S*~(RQ»SQRT(RAn))/(?.C*AQ) 

X«XSAT»ELO*S 


_ 


5J> 


non 


t r v'  i 

i 


LATLON 


DATF  * 76261 


16/02/47 


I 


7*/SATfcN0*S 
PRINT  20 

20  FORMAT I / / 2X ♦ *L  AT LON  : CL OHGT , RASQ t AQ, BQ,CQ tR AO t S, X * Y» l • ) 
PRINT  ?l,CinHr,T,BASC,AC.BCtCGtBACtS.XtY»  l 

21  FOPMATI /2X, 10F11.4I 


CONVERT  TC  EARTH  COORDINATES 

I 

XL AT= AT «N( 7/SORT (X+*24Y**2I) 
XLON  = AT  AN2  C > » X ) 
XLATaXLAT/ROPDG 

XLCN=XLCN/RDFOG  

GO  TO  40 
32  NF  R R*  2 
40  PFTURN 
END 


*, 


PPECFS 


0 A Tf:  = 7f»26l 


1 6/ 02 /4  7 


7FI  21 

SU«*>nUTTNP  rPFf.FS  ISAMTIM,SPAX1,SPAX2,SPAX3 ) 

rm-MON/SC  AM”  71  SCAN  ♦NLMSFN  *NOPC  LN  fTCT  L 1 Nf  QPOL  T H,  RADI  IN 

l » PICL IN, TOT  PI F» CFGPLF»PADELE,PICFLE ,FF,PI TCH  , YAW ,ROl l , SKFW ,R0TM1 I 
2.RDTR13  t»HTP2l  ,RCTM23»PCT*3I  »R0TP33 . RF ACT « RC AS  IN »S C. CO, P 0 IR . PR AT 
if  (prat.fo.o.oi  r,n  to  10 

PTOT=S APT IM*  PR  AT 
?A  = $lNlPTCm 

r.A=r.nsi  ptot  » 

XI  =S(JRT  C 1 -0/ ( 1.0M$PAX17SPAX3I**2  I ) 

Y 1=0.0 

n = -(SRAXl*Xl/SPAX3l 

X2=S  PAX2*7l 

Y2=SPAX  3*Xl-SPAXl*Zl 

72=-SPAX2*Xl 

Xl=CD*XU$0*X2 

Yl=C0*Yl»S0*Y2 

Zl=CD*7 1*SC*Z2 

PRINT  20»Xl,YI,Zl,X2.Y2tZ2 

20  FCRMATI//2X, ‘PRFCES:  XI  ,Y  l ,Z It X2, Y2, Z2 •» /2X,6F 15. 71 
$PAXl=CA*SPAXi»SA*Xl 
SPAX2=CA*SPAX2+SA*YI 
SPAX3=CA*$PAX3*$A*Z1 
10  RETURN 
FNO 


G IPVFL  71 


L INTLF 


OATF 


7ft  2ft  1 


1 ft/02  /ft  7 


| 


MjnonijTINF  LTNriF(XF,Yf:,/F,XllN,XFLF) 

COMMON/ SATVFC/X SAT  , YSAT  ,7  SAT  ,HF  IOHT  , XVECl  . XVEC2  , X VEC 3,  THC T A , 

* AROTll  , AR0T12,  APFT13,  AROT21,  AROT2?»  ARnT23,APOT  31,  AP.OT17,  APOT3  3 
COMMON/ SC AN»/ 1 SCAN  , NllMSFN  ,NOPCLN ,TOT L l N, OEGL l N«  P ADL  IN 

*,PCI N , TnTELc,rFGFLFfo AOELF, PIC FLE.FF, PITCH, YAW, ROLL, SKFW,RQTM1 l, 
♦POTMn,PnTP21,RCTP2  3,FCTM3l  , ROT R33 , R FAC T , RO AS l N, S 0, CC, PO IR , PP AT 
Y 1 = XC*«X  S A T 
Y2  = Yr--YSAT 
Y3=ZE-2SAT 

Y F A C T = 1 .0/ SORT  I Yl*«2+Y?**2*Y2**2> 
yi=yi*yfact 

Y2-Y2'*  YFACT 

Y3*Y3*YF ACT  

XI =APOT l l ^Y1 ♦APCTl?^ Y2+ARCT13*Y3 
X2*AP0T2  1*Y  l + ARnT22*Y2«,AR0T23*Y3 
X3=AP0T?l vY1*APOT32*Y2*ARCT33*Y3 
UMV=ATAN2( X3,S0RT(RFACT-X3**2J l-RCASIN 
XL IN=PCLN-UMV/RAOLIN 
SI  IN*SINUJPV» 

CL  I N=CO  S ( LM  V ) 

U=R0TM1 1*CLIN*R0TM13*SLIN 
V=RnTM2l*CLIN*PCTK23*SLIN 
PRINT  5,U,V,UMV 

5 F0PMATI/2X, ‘LINELE:  U, V,U*Y* ♦ /2X , 3E20. 8) 

IF( V.PO.O.O)  GC  TO  6 

UV=ATAN2(V,U) 

GO  TO  7 

6 UV=0.0 

7 UMV=UV-ATAN2(X2,X1) 

XELF=PICFLE*UPV/PAOELE 
PRINT  10 

10  F0RMAT(//2X,'LINELE:  Y F AC T , Y 1, Y2 , Y3, Xl , X2 , X3 ,UMV,U, V» ) 

PRINT  11, YFACT,Y1,Y2,Y3,XI  ,X2 ♦ X3 ,UPV ,U, V 
. 11  FORMAT! /2X,7Fll.ftt/2X,3E20.8) 

RETURN 

END 


i 


t 


* 


| 


ono  non  nfinnnn^nrt'i'^ 


LEVFL  ?\ 


VCLALO 


D AT  F •=  70?. ft l 


16/02/47 


SUPROUT 1NP  VCLALHf  ATFR AC , XL  AT , XLCN, XE , YE , Z E t NT0 R ) 

RE  ATE  Rn()  = I350,6)SATPAK,FTK/NV 

X F » Y C » 7 F - COMFCNT-NTS  OF  EARTH  COORDINATES  VECTOR 
YLAT  = X L AT  CCNVEOTFD  TO  P A C IANS 
YLCN  = XLTN  CONVERTED  TO  RADIANS 
SINLAT  * SIN  ( YLAT  ) 

COSIAT  = Cf'  ( YLAT  ) 

S I ML ON  = SIN  ( YLQN  » 
enSLON  = COS  ( Y LON  > 

C.  i . MOTTCRSHEAp/CSC 

CCMMCN/GDATA/PIfPCPCG.R, A, B . AB , ASCt BSQ, ATMHGT.GR ACON, EMEGA , SOLSI D 
* « SHA  * I YR  » I HP 

COMMON/ SA  TVFC /XS AT  »YSAT»ZSAT, HEIGHT*  XVEC1 *XVFC2»XVEC3*  THFT  A , 
*ARDTU  , AR0T12,A90T13,AP0T21,  AROT22,  AROT23,  AROT3 1 , AROT32 , AROT 33 
YLAT*XLAT*<RCPDG 
YLCN*XLCN*RCPCG  . 

SINLAT*SIN(  YLAT)  

COSLAT=COS( YLAT) 

SINLCN=SIN(YLON) 

COSLCN*COS ( YLON) 

X«COSLAT*COSLON 
Y»COSLAT*SIALCA 
2*SI NLAT 

CHECK  IF  POINT  IS  CUT  OF  SATELLITE  VIEW  AND  IF  SO  SET  ERROR  FLAG 

IFI  ( X*XVEC1*Y*XVEC2*  Z*XV,EC3I  .IT. THETA)  NERR«3 

; , 5 

ADJUST  FOR  CBIATENESS  OF  EARTH  SPHERE  /yND.CLCUD  HEIGHT 
TANLAT*!  SIHAT/CCSIAT  )**2 

. RR*SORT( ( L.O*TANLAT)/(BSO+ASQ*TANLAT) I*AB*ATPHGT*ATFPAC 
XE«PP*X 
YE*RR*Y 

ZE»RR*Z  _ . _ ■ 

RETURN  ' I 

END 


I 


J 
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FUNCTION  VT  IMF  VICTIM,  XI  IN  » XFLC  ) 

CrM«nN/NAVSl6/IMV,N/lVN,L£Nf:N,NIT*MTf  NORrt.NHAV,  FL»  F.P,  FT,  SP  INC  3 I , 
l RASf.FN,  OFfL  IN,SPINRA,TMPSCL  ,GT  I Ml  16 » »BFT  A (16)  ,«DOT  I 16 ) , MG  AM  U 6 ) 
r.r^MON/SCANR/  IS  CAN  f NUMSEN  *NnPCLNt  TOTL  I N*OFGL  I N » R AOL  I N 

• fPICLIN,  TnTFLEfOFGFLF  LF  ,PICFLF»EF»  PITCH,Y  AW,«>Oll,SKFW  ,ROTMl  1, 

*ROTMi3fpnTR2ltPOTM2.1fROTM31lRnTR?3*RFACTtROASlN«50fFD*POIR  « PR  AT 
TLlN»XLIN*-0«5 
PAR  L IN= I ILIN-ll /NtMSFN 
PAOFLF= (XElE-l.O)*EF 
FRAMFT=TMP^Cl*(FAPLIN*PARFLr) 

VT  IMF=P ICTIRfFRAMFT 
PRINT  10,PAPL INfPARFLE, FRARET.VTIMF 
10  F0DMAT(//2x,*VTIMEs  PARLIN,PARFLE,FRAMET,VTIMF' f/2X, 6F20.8) 

RETURN 

FNO 
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Suppoutinf  rfsitu 
CrMMPN/r  YSf.C**/  |T«fM 

rpvMtiN/XL  ANO/NLANO,LOAV,I  fnnri  32)  ,PT  I MR  ( 72  ) , XI  I N|32  I , XEL  F (3?  ) , 
l*lATn2l,XLrNI3?l,CLlN<3  2 J,OELTA<32),TIMPLl  32| 
CMMMON/NAVHN/IKAVtNAVK  , L A ADA  , MT  , F I T , NORR , NO  AY , FL  , FP , FT , $P  INI  1 1, 
1 R ASC^O'CL IN.SPINRA,TMPSCL,GTIM| 16) ,RPTA< 16) ,ROOT( 16) ,NGAMI16) 
n I MF NS  I r N NUN  1 16 ) » AVG ( 16 ) , VAR  1 16  ) 

G AM A=0. 0 
GOPT-O.O 
ATPR.0.0 
PO  10  L«1.W 

NUMILWO 

AVGID-0.0  ..  ..  ... 

VARID-0.0 
10  CONTINOF 
LMAX-1 

PPINT  17,  NC AY 

17  FORMAT! / /20X  , 68HL  ANCM ARK  RFSICUALS  AT  ZFRO  GAMMA  SHIFT  FOR  OA 

1TF  ,16) 

PRINT  18 

18  FORMAT  1 2 X , 1 10H  IMK  RETURN  TIFF  I HRS ) LATITUDE  LONGITUDE  MEA 
IS.LINF  CALC. LINE  LINE  ERR  MEAS. PIXEL  CALC. PIXEL  DPIXEL) 

SUMSQ=0.0 
7L  SQ=*0, 0 
PRINT  17 

37  FHRMATIUX, 'LANDMARK  FESICUALS  AT  ZERO  GAMMA  SHIFT', /2X,  ' LMK  GMT 
1 LATITUDE  LONGITUOF  CALC. LINE,  ERROR  CALC. PI XFL, SHIFT  RET • ) 
OP  50  I * L , NL  AN  C 
PICTIM-PTIMFI  I ) 

TIMFH  I l-VTIME  I PICTI  M , XL  INI  I ),XELE(  I ) ) 

NFPP.  *0 
MSER  * 10 

CALL  SATPAR|PICTIM,YL  IN  , YE.LE  , XL  AT  1 1)  , XLONII ) , 2,  NFRR,  GAMA  , 

1 GOOT,ATFR) 

IFINEWR.GT.O)  MSER*20 
DELTAII)«YELF-XELE(  1 1 
DLINID-YLIN-XLINII  ) 

PRINT  23,l,MSER,PTIME(I),XLATI I), XL ON II),XLINt!),YLIN,OLINfI ), 

1 XFLEII ) ,YELE ,DELT  A ( I ) 

23  FORMAT! 2X,I6,2X, 16,  3F 1 2 . 5 ,6F 1 2 . 3 ) 

ICOOENaMOdl ICOCEI 1)7100,10) 

IF! ICOOEN.NF.O.ANO.ICCCFN.NE.l  ) GC  TO  60 
SUMSQ»SUMSQ*DELTAII )*OcLTAI  I ) 

ZLSC*  ZLSC*CL  INI  I )*CL  INI  I ) 

LMK  «MOP 1 1 CODE  I I), 100) 

IFILMK.LT.O.OR.LMK.GT.32)  LMK» 1 
AVGILMK) «AVGILFK) *CL  IN  (I ) 

VAPILMK)  «VARILMK)*0LIM  I ) **2 
IFILMK.GT.LMAX)  LMAX*LFK 
NUM(LMK)aNUNILMK)*l 

60  PRINT  67,  I ,PT IMEI I ) , XL  A T I l ) , XLCN ( I ) ,YLIN«Ol INI  I ) , YELF , 

1 OrLTAII ),NFRR 

6 7 F0RMATII6,F7.2,2F10.3  ,6F20.6,I6) 

50  CONTINUE 
, FNIJM  *L  ANON 

;ms«sop.ti  zlsq/fnum) 

RMSaSQRT ISUMSQ/FNUM) 

PRINT  5 7 »RF  J , ZFS 

i.f. 
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S7  FORMATION  OPS  PIXEL  FPR*  ,FH.2,16H  RMS  LIMP  FRP»  ,FR.2I 
POINT  M * 7 M S 

61  FOOWAH*  summary  OF  SPIN  ATTITUDE  fit:  OVERALL  RMS  LINE  ERROR » • 
l »Fe.2»nxi 

no  so  lmk*i,lmax 

IF  IMIM  IRK)  .EC.O)  GO  TO  SO 
FNL«=NUM  LVK) 

QLAV=AVG<LMK) /FNUM 
VSC=VARILRK)/FKUM-nL«V*'»2 
VSQ  = SQR  T ( VSO I 

PRINT  71,  LMK,  CLAV,VSQ,MJM(LKK  I 

71  FORMAT ( • LANDMARK* , 16, ' MEAN  LINE  ERROP F8.2 F8. 2, ', ON, 

1 16, *,I MAGFS.' ,5X1 
90  CONTINUE 
EL  *Z  M S 
RETURN 
FND 


* 
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SIIPROIIT  INF  GAMCAL 

C TMvr  N/  XL  ANH/NL  ANC  #L  PAY » 1C  CPF  ( ? 2 » , PT  I M*  ( 32 ) , XL  I N(  32) » XELF ( 32  I » 

1 XL  AT (32) ,XL0N(32)  ,OLIN(32) ,DF L * I 32 ) , T I *EL I 32  ) 

COMMPN/NAVSLN/I  KAV«NAVN*LANCN*NI  T * M f TfMnRfUNfHYt  FL,EP,CT,$PIN(3), 

1 PASCFN,OFCLIN» SPINPA,TMPSCL,GT 1M(16),BETA( 16), BOOT ( 16),NGAM(  161 
CALL  GCOOER 
PFTUPN 


END 
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SlIHRnilT  INC  GSM1  FT  ( MJVf  AM  ) 

COMMON/ XL  ANn/NLAND,LDA Y,  ICCDE (32)  ,PT!MF(  32)  ,XUN<32)  ,XFLFI  32) « 

1 XL  AT  (32  ).XLON(32),CL  1M  3 2 ),DFir  ( 32  ) , T1  Mf?L  ( 32) 

CnMMnN/NAVSlN/TNAV,N/VN  ,l/NCN,NtT, PIT, NORR,NOAY, CL  » FP  » FT  , $P  IN(  3 ) » 
1 PASr'  N.Oer.LlN.SPINRA.TMPSCL.OTIMl  16 ) , BET  I ( 1 6 ) , BETD(  16  ) , NUMG  ( 16  ) 
INAV=1 
NUMG  AM a 0 
! ST  AR  T*  1 

IF  ( ISTART.GT.MAAD)  GC  TO  6 

SUMG =0.0 

SUMGT*0.0 

SUM  T =0.0 

SUMTSQaO.O 

NUM  = 0 

N*0  , i 

no  s i= i start ,nland 
IFII.A4F.IST/RT)  GO  TO  13 

PT!M=PTl*c( i | 

TCODEN^PODI I CODE ( I ) / 100 , 10 ) 

icodfg=pco(  rccoem/iooot  ion) 

!F( ICOOFN.NF.O.AND. ICCDEN.NB.2 ) GO  TO  2 
AUM  = NUM«-1 
GA  MMA=DE IE ( I ) 

SAMT  IM=TI MFL ( I ) 

sumg=sumg*gamma 
SUMGT=$»JMG7*GAMMA*SAPTIP 
SUMT=SUMT*SAMTI  M 
SUMTSQ  = SUMTSC+S APT  I P** 2 
PRINT  40 

FORMAT (2X, 'NUM,  GAMMA, SAMT IP, SUMG , SUMGT » SUMT , SUMTSO ' ) 

PRINT  41  ,NUP,GAPPA,SAPT  IM  , SUMG, SUPGT, SUMT, SUMTSO 
F O P. M A T ( 2X,I8,3E20.4,/2X»3E20.4) 

I F ( I .FQ.Nl ANC)  GO  TO  20 
NCOOF*MOO( I CODE  I 1+11/1000,100) 

IF(NCODE.FQ.ICOOEG)  GO  TO  5 
N«I 

IF(NUM.FO.O)  GO  TO  10  

nu«gam*numgamm 

GTTM(NUMGAM)=PTIP 
NUMG (NUMG AM ) «NUP 
I F ( NUM. GT  • 1 ) GO  T (f  4 
BF  T I ( NUPGAP ) *GAPMA 
BFTOCNHMGAM  |«0.0 
GO  TO  10 
XNU*»NUM 

OFMOMa  XNUM ♦ SUMT  SO-SL’M 7**2 

BE  T I (NUMG AM ) « (SUPTSQ+  SUMG- SUMT  +SUMGT ) /DENOM 
RETO(NUMGAM')  «|  XM.M*SUPGT-$UMG*SUPT)/DENCM 
PPINT  50 

FORMAT ( 2X  ,»GSHI FT : XNUM , DENG* , eET l , BETD * ) 

PRINT  XNUM  ^ DENOM, BE  Til NUMGAM) ,BETO( NUMGAP ) 

FORMAT ( 2X» I8.3F15.7I 
GO  TO  10 
CONTINUE 
l $ T AP  T »N+ I 
co  T0  l 
. PFTUFN 
END 
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subroutine  Gf.nnFR 

(TMMPN/SYSCOM/  ITKfNL 

r rf  MON/XL ANO/NLAAn.inAY  , TC0CEI32J.PT  IMFI32),XLINI32)  »XFLH  32), 

1 XL  AT  ( 32  ) . XLONI  22 ) ,UL I M 32 ) .OELF ( 22 ) , T I MEL ( 32  ) 

CFMMC  N/NAVSLN/  I NAV.  NAVA.  LA  NON,  N l T , M l T ,NORR ,NOA Y,  FL»EP»ET«SP!N(3)» 
l R A SCEN  . DEC  L IN  , $P  I NR  A ,TN  PSCl  ,GT IMC 16 ) , BET  I ( 1 6 ) , PET D<  16  » , NUMGI 16  ) 
riMTNSirN  RF I T ( 32  1 » POLO  I 3 2 » » RNEW  ( 3 2 1 , 1 MI  D ( 32  I tLMK  I 01,32  I 1 1 HD  ( 16 1 

0 IMF  NS  I CNR  MS  1 16  I , JCCCE ( 32  ) , LAST < 16  ) 

INTFGEP*A  I T I ME ♦ I L ALO 

NLP  » 20 
PP=0.0 

NPAGEsl+NLANO/NLP  ' 

1 MG*  1 , _ __ j 

TIMAG=PTJMF(ll *0.01 

00  10  1*1 . NL  AND 

lEfPTIMEm.LT.TIMAGI  GO  TO  6 

LASTCIMG)*I-1 

IMG*TMG*l 

T IMAG=PT Imc| j) f0.oi  ' 

6 I M 1 0 U)  * I MG 

LMKIOU  l = MOCncnCE*  I),  100)  - 

RFlTm*0.0 

JCO0E1 I I * ICCOE ( I I 

pOLom*o.o 

RNFVMI  1=0.0 
IW0(  I )*NLANDM 
10  CONTINUE 
EP*O.0 
IGP=  1 

MAXI M*I MIO(NLANC) 

LAST! MAXIM  )*NLAND 
IW0(1)=MAXIM 
60  NRGN*1 
. 00  90  N* 1 . NL ANO 

T M* I W0( N I 

IFIIM.GT. MAXIM)  IM.MAXIM 

_ IF  (IM.GT.O)  GO  TO  70 

PRINT  20 

20  FORMAT ( // * NEGAT I VE  ENTRY  TO  GAMMA  CODE  ENCOUNTERED*) 

70  N«=NO*LAST  UM) 

00  60  I *NBGN  fNEND 
ICOOFC I )*JCOOE (1 )+lOOO*N 
80  CONTINUE 
NBGN*NEN0*1 

IFINBGN.GT.NL ANC)  GO  TO  95 
SO  CONTINUE 

95  CALL  GSHIFT  INUMGAM) 

NUMG ( 32 ) *NUMG A M 
NTOT  *0 

OP  99  I GP* 1 .NUMGAM 
NTPT*NTOTfNLMG( IGF) 

RMS ( IGP  )«0.0 
99  CCNTINUE 
KHGN*l 
KEN0*0 
T0^S0«0.0 
PRINT  701 

701  FORMAT  I 16X  » A 1H**  LANDMARK  PIXEL  SHIFT  COMPENSATION  **  ,17X) 
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PRINT  707 

707  FORMAT  ( ' TC  CHAT  PCI  THE  E7ST-WPST  AL  IGMrNT  nF  THIS  S^T  OF  M/VSTF& 

1 IMAGES*  • »/2X,  »FXAM!Mr  7MF  MF  A SLpF  0 LANDMARK  PIXEL  SHIFTS*  AND  OFF 

2 ! Nc  IMAGF  GROUPS', /2X,  'FOR  A TP  I AL  LINEAR  FIT  OF  SHIFT  VS.  TIME.  T 
3HC  PAR  AMFTFPS  ANO  F 1 N AL ' ♦ l X , ' ( NEW ) PFSIDUALS  FOP  THF  CURPFNT  GROUP 
4 I N G AP.f  LISTFO  RELOW:  •) 

PRINT  702 

702  FORMAT!  • MASTER  IMAGF  LANDMARK  SCAN  PIXEL  SHIFTS  FIN 
I AL  PIXEL  ERROR  •) 

P’INT  703 

703  FORMAT!  • IC  GMT  GROUP  NO.  10  T I MF  MEASUREO  F ITT EO  • 

I,'  NEW  CLO  •) 

N*  I 

18  LFN0=0 
K P AGF  = 0 

19  LBGN  = LCND*-1 
LEND=LEN04NLP 
KPAGE=KPAGE4l 

IF!LENO.GT.NLANC)  LENC=NLAND 
IGP=0 

22  DO  50  I = L BGN  *L  ENC  , ' l 

ITIMF*ILALO(PTIME!IH 
IGOLl)=  I GP 

IGPsMOO! I CODE! 1 1/1000, LOO) 

BETA  = DE  T I ( I GP ) 

RETOOT  sRETC(IGP) 

30  ROLD  ( I ) sRNFWt I ) 

GCALC=BETA40ETD0T*TIMEL<  I > 

RF IT ( I) =GC ALC 
01 ff  =DELE ! I J-GCALC 
RNEW!  11=0  IFF 

IC0nFN  = wnO( ICCOE!  I)/ 100, 10) 

IF ! ICOOEN.NF.O.ANO. ICD0FN.NE.2)  GO  TO  40 
. CSOniFF+"-2 

RMS! IGP) =RMS( IGFJ40SC 
T0TS0=T0TSC4DSQ 

40  PRINT  704,  IMID!  I),  ITIMF,  IGP,I,LMKID!I),TIMELm  , 

IDELFI I ) ,BFI T! I ) ,RNFW( I ) ,ROLD!I) 

704  FORMAT! 14, 18, 14, 16, 14, 5F 18. 4) 

50  CONTINUE 

T0TS0*S0RT!  T(ST  SO/FLO  AT  ( NTCT  ) ) 

PRINT  708,  TOT  SO , EP 

708  FORMAT! 32 X, • FINAL  RMS  PIXEL  EPRCP-' ,2E20.4) 

FP=TOTSQ 

fTN*$CRT! (EL#+24EP**2)/2.0) 

PRINT  709,  ETN,ET 

709  FORMAT! 12X, 'CORRESPONDING  OVERALL  NAVIGATION  ACCURACY  *,2E20.4) 
FT«5TN 

PRINT  219 
219  FPRMAT174X) 

PRINT  77 

77  FORMAT! /5X» 'GAMMA  SHIFT  CALCULATION*  ) 

PRINT  230 

230  F0PMAr|«4H  GROUP  SIZE  PASETIME  BETA  BETADOT  RMS  I 
00  250  N*  I , N(JMG  AM 

PMSIM  «SCRT!RMS!N)/FLCAT  INUMG(N)  )) 

PRINT  240,  NtNUMG(N) .GTIMINI ,BETIIN) ,BETO(N) ,RMS!N) 

240  FOP**AT  (216, 4E20.4) 
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250  CONTINUE 
700  Kf  TURN 
FNO 


APPENDIX  E 


Landmark  Scanning  Program 

by 

Neil  R.  Guard 

Rufus  E.  Bruce 


Sandra  K.  Weaver 


6-I5.*2b:ss  (,0) 


COdE(I)  0006R7J  o AT*  ( 0 1 OOi(l3o:  BLANK  C<jMM0N(2)  OOOOOO 

BENCES  (BLOCK.  NAME) 

P 

s 

i 


*ment  (Block,  type,  relative  location,  name  i 


1 

1 1 1 G 

ouni 

000 1 73 

l R R G 

000  1 

0002 | R 

1 72G 

OOqI 

0(3(3  2 2 2 

200G 

j( 

2 I 7 g 

0001 

000  3 | 2 

236G 

onoi 

000B5A 

2776 

OUoi 

(J(J(3  R 6 3 

3 (J  3 G 

> 

331  G 

0001 

ooossn 

3r3G 

0001 

000566 

35RG 

OOoi 

000636 

3 7 7 G 

3 

7S0L 

0001 

000Rr*3 

7“7SL 

0001 

0006  1 3 

810L 

000  1 

0006  1 7 

850L 

3 

> 

sjOf 

0000 

00B0A2 

920F 

oooo 

00R03R 

V30F 

OOoO 

00R036 

9 R OF 

960F 

oooo 

00** 0 j 2 

970F 

00(30 

I 

003773 

1 

OOOO 

I 

003766 

I AFC 

♦ 

IDMP 

OOQO 

l 

Q03/,  36 

tELE 

oooo 

I 

003771 

1 EL  1 

OOoO 

I 

0Q3772 

IEL2 

7 

IPAr 

oooo 

I 

003760 

ISKIo 

oooo 

1 

OOROOI 

ISLP 

OOoO 

1 

00R002 

ISLPI 

r 

ITOS 

0000 

I 

OOSfrS 

I VAL 

oooo 

I 

0037RR 

1 RO  1 

oooo 

I 

003752 

1 RO  2 

> 

I 

1 X I T 

OOOO 

I 

00377R 

J 

oooo 

I 

003765 

JOUMp 

00(30 

I 

003775 

JX 

♦ 

line 

OOOO 

I 

001  3rR 

nata 

oooo 

1 

OOROOO 

NOX 

0000 

I 

00R003 

nline 

» 

NpOrdI 

oooo 

I 

0037a! 

nnoros 

oooo 

R 

0002 | R 

T 

OOoO 

R 

002R7R 

X 

c •••< 

c 


c • • • • 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c •••• 


this  program  rfads  a sectored  version  of  an  sms  data  tape 

AND  PROCESSES  A USrR  CHOSEN  AREA  BY  USING  A PATTERN  RECOG- 
NITION SCHEME  TO  LOCAtp  A PREDETERMINED  LANDMARK  FOR  USE 
IN  CALCULATING  REGISTRATION  TRANSFORMATION  PARAMETERS 


initial  STORAGE 
1TOSU3R) 

K ( 6 > 

T(60.10l 
NATA(60»I0)  . 
X(AO.IO) 
LINE(IP) 

I ELF-  ( 10) 
ICRTELUO)  . 
I WO  1(6) 
l*D2(*) 


data  return  from  flotap 
status  return  from  FLOTAP 
COASTLINE  OUTPUT  ArraT 

forking  storagc  array 
slope  contrast  array 
CRITICAL  line  NO,  STORAGE 
CRITICAL  ELEMENT  no,  STORAGE 
COASTAL  element  STORAGE  no, 

character  computation  storage 
CHARACTER  COMPUTATION  storagc 


c • ••• 

c 


DIMENSION  ITOSi t SN I ,K  ( 6 ) ,T(60, |0| ,NATA(40,I0) « X I Ao . | 0) , LINE  I 10) , 

IELE ( 10) . ICRTEl ( AO) , I WO  I (A  ) , 1*02 ( 6)  ■ 

INITIAL  PARAMETERS  73 


% 


C I SK  I P - No.  OF  Llrirs  ON  TAPf  to  SKIP 

C MAOKOS  - No.  OF  WORDS  PF  R L>NE 

C NLINES  - No,  OF  LI  NFS  IN  ARRAY  ' 1 

C IPHP  - No,  OF  WORDS  TO  SKlp  AT  BEGINNING  OF, EACH  LINE 

C JDMP  - No,  OF  WOhoS  LEFT  To  SKIP  AT  END  OF  EACH  LINE 

C IAFC  - For  SAvINg  NO.  OF  ABNORMAL  FRAME  COUNTS 

c ipar  - for  saving  no.  of  parity  errors 

C •••• 

I S K I P ■ 200  1 

NVOROS-  1 0 , j 

NttOROI-NWOROS-l 
N L I N E S ■ 60 
I DMP«7 I 

JDUMP® I0MP*NW0rdI“1 
I AFC"  0 
IPAR»  0 

ILJNF-  ISKIP  ♦ NLINFS 
I El.  1 * I OMP* 6 
IEL2-  I EL  1 ♦ NwioROMG 

c • skip  a page  on  output 

WRITE  16,9001 
900  FORMAT! I H 1 ) 

c • establish  tape  pitmen  program 

CALL  FLOTAP!  I I , .VlSlBLEOATAl  « I 
CALL  F L 0 T A P C 6 )' 

CALL  FLDTAPiS! 

CALL  FLOTAP! I I 

C • SKIP  ISKIP  LINES  oN  TAPE 
CALL  FL0TAPI9,  tsK|P| 

C • LOOP  TO  REAO  PATa  AND  STORE  IT 
WRITE  (6,970) 

•970  FORMAT  UX,*RAp  DATA  ARRAY  AS  READ  PROM  TAPE**/) 

WRITF  ( 6 ,-950  ) MKIP,ILInE*!ELI  ,|EL2 
950  FORMAT  C I X , • L I NF * . I A , * TO*, 16,*,  ELEMENT *,  I 6 , * T0*,|6,//) 
oo  ion  i*i  ,nl Infs 

C • CLEAR  FLOTAP  STATUS  RETURN  ARRAY 
DO  200  J«  I ,6 
K(J)«0  e 

200  CONTINUE  * - ...  - 

C • PICK  UP  DATA  LINE 

CALL  FlDTAP(7.13B,IT0S,k) 

C • CHECK  tape  read  status 
IF  IKI3J.NE.0)  STOP 
IF  ( K ( 9 ) ■ NE , 0 ) STOP 
IF  (K(S).NE.O)  STOP 
IF  (K(I).EQ.O)  STOP 
IF  JKI2I.NE.0)  |AFC«IArc*l 
IF  <r<6)*.NE.O)  j P AR • l P AR ♦ l 

IF  <K(2).nE.O.anO.K(6),nE*O.ANO.K< I I.EO.l ) STOP 
C • OUTPUT  DATA  LINE  READ 

WHITT  16,9301  ( |TOS( JXj , jX - | OMP , jOuMP ) 

930  FORMAT  ( | X ,9A6) 

C • STORE  desired  elements 

DO  900  Jr  I ,N*0rdS 

JX-|DHP*J-I  74 

NATAI I ,J)*lTOS(jX) 

900  CONTINUE 


5 


* V 


p 


loo  CONTINUE 

WRJTF  16,900) 

C • LOOP  THROUGH  NLlMF  LIMES  TO  process  OAT* 

DO  son  1 = 1 ,NL  Infs  (U) 

I X I T ■ 0 

C • STORE  SINGLE  ELEmfMT  VALUES 

00  600  J=  I , Inform  i)f\ 

1*01  ( I )=FLD<nO,*,NATA<  I , j)  ) 

1 *0  1 (2)>FLD(06,*,NATA( I , J)  ) 

1*0113)  =FLO(  12,  *,  NATAll.J)) 
Ittl)l(9)=FLD(|8,6,NATA(I,J>) 

1*01  (S)eFLOI  29  , , N A T A ( J , Jl  ) 

16  01  ( 6 ) *FLO ( 30  , f ,NATA I ) ,J)  ) 

I 6 0 2 ( l )«FLO( |0,6,NATA(  j , J ) I 
IWD2(2)«FLD(29,/,,|JAtA(  ( ,J>  ) 
IWD2(3)=FLO(30,4»NATA( | ,J)I 
IW02t9)aFL0(00,*,NATA(  I ,J*l  ) ) 
I*D2(5)=FL0(06t*,NATA( j , J*  1 ) I 
I*D2(6)*FLD(|2,'a,NATA(I.JM>> 

C • FIND  SLOPE  ACROSS  FOUR  ELEMENTS 
00  700  I X ■ I ,6 
NO X ■ 6 • ( f X- I ) 

I5LP*  ABS I 1*01 ( ) X ) -1*02 ( I X ) ) 

C • SET  OCEAN  VALUE  f r SLOPE  ,lT.  3 
IP  (ISLP.LT.3I  SO  TO  ?2S 

C • SET  CLOUD  VALUE  IF  SLOPC  .ST.  9 .Or.  ELEMENT  .GT. 

IF ( ISLP.GT.N.Or. I 70) C I x ) .6T.22)  GO  TO  750 
C • SET  COASTAL  VALUp 

FLD(MDX,6»T< I ,J) l«29 

IX|T«IX!T*I 

GO  TO  775 

7S0  FLO ( NOx , 6 , T I I ,j) )«8 
GO  TO  775 

725  FLOINOx ,6,T( I , J ) ) >5 
C • STORE  SLOPE  VALU£  FOR  OUTPUT 
775  FLOINOX, A, XI  I ,j)  )=  |SLP  ♦ 98 

IF  (IXIT.E0.1I  i CRTEL  < I ) ■*• J4  I x + 2 
700  CONTINUES  , . _ 

600  CONTINUE  *1  V \ * t 

500  CONTINUE  Cf 

c • output  coastal  outline  array  ' 

WRITE  (6,990) 

990  FORMAT  ( IX. ’PREDICTED  COASTAL  OUTLINE  ARRAY*, /) 
*R|TE  (6,950)  ISKIP.1LINE.IELI.IEL2 
00  110  I ■ I .nlines 

WRITE  (6,9301  ( t ( 1 , J ) , Ja  | , NvvORo  I I 
I 1 0 CONTINUE 

WRITE* 16,9001 

c • output  slopes  array 

WRITE  (6,9601 

9*0  FORMAT  (I  X , *COM»*,tTEO  SLopF  VALUE**,/) 

WRITE  (6,9501  |SK  IP,  ILluE • (ELI , IFL2 
00  170  I ■ I , NL  I NF 5 

WRITE  16,9301  ( X ( I * J ) , J= I , NWORO I ) 

I 20  CONTINUE 

WRITE  (6,900) 

WRITE  16,9101  ltCRTFL(|t,|a|,60)  ic 


C*l  VI'1 

Cl 


(,  . \ 


9)0  FORMAT 

• calculate 


I IX,  'CRITICAL  ElFMENTS*  ,/,6(|X.I5)  » 

L I Nf  Anr,  ELEMENT  NUMBCRb  FOR  UP  TO  TEN  COaSTAL 


C PROJECTION  point?;. 

ISLPI-  ICWTEI  ( U- ICRTFli  3>  ' 

J-0 

NL  1 NF I »NL I NFS-2 

l>0  800  I »2  » NL  1 NF  J 

15LPZ=  ICRTELI  1 l-ICRTELI  1*2) 

I V AL " I SLP  I • I 5LP2 
IF  (IVAL.6T.0I  00  TO  6in 
IFU.GF.IO)  GO  jo  850 
J»  J*  1 

L I NE  I J I * I ♦ 1 
1ELE(  J)a|CRTEL(  i + n 
8 | 0 I S L P 1 * I SLP2 
800  CONTINUE 
85O  JX»J 

c • output  critical  ljne  and  element  numbers 

WHITE  16,900) 

WRITE  16,920)  lLlNE(l),|ELE(l),I«l,JX> 

920  FORMAT  ( IX, 'COASTAL  PROTRUSION  POINTS  PREDICTED  *x 
1 IOC/, 91X, -'LINE  *,12,*  ELEMENT  M2)) 

STOP 

END 
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